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Abstract 

We extend the existing theory of approximation orders provided by shift-invariant subspaces 
of L 2 to the setting of Sobolev spaces, provide treatment of L 2 cases that have not been covered 
before, and apply our results to determine approximation order of solutions to a refinement 
equation with a higher-dimensional solution space. 
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1 Introduction 

1.1 General 

We are interested in this paper in the approximation order of shift-invariant (SI) spaces of func¬ 
tions defined on the Euclidean space IR'^, d > 1. Such spaces play an important role in several 
areas of Real Analysis, including Spline Approximation, Wavelets, Subdivision Algorithms, Uni¬ 
form Sampling and Gabor Systems. It is not surprising, thus, that the theory of approximation 
and representation from SI spaces received significant attention and enjoyed rapid development in 
the last 10-15 years. The determination and understanding of the approximation orders of these 
spaces is among the main pillars of this theory. 

As the title of this article indicates, we restrict our attention to approximation in Sobolev 
spaces: given s £ IR, we denote by VU2(1R'^) the Sobolev space of smoothness s, i.e., the space 
of all tempered distributions / whose Fourier transform is locally in L2(1R'^) and satisfies 

+ I ■ I) 

(Here and elsewhere, | • | is the Euclidean distance in IR'^.) A closed subspace S C Wl(IR'^) is shift- 
invariant if it is invariant under all shifts, i.e., integer translations, or more generally, scaled 
integer translations: given a fixed h > 0, 

for every a £ hTL^ and every / £ lUI(IR'^), f € S /(• -|- a) £ 5. 

When necessary, one identifies the underlying parameter h by referring to S as /i-shift-invariant, 
and/or by denoting the SI space as S^. Also, sometimes, in order to emphasize the ambient 
space VUI(IR'^) we write instead of simply S. The smallest SI space that contains a given 

$ C lUKlR'^) is denoted by 

5$:=5$(W|), 
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or, in complete detail, S^iW^), and we refer then to as a generating set of 5'$. The basic 
objective of shift-invariant space theory is to understand the properties of SI spaces in terms of 
properties of their generating sets. In this regard we recall that an SI space generated by a singleton 
‘h = {i^} is known as Principal Shift-Invariant (PSI), while the one generated by a finite ‘h is 
referred to as Finitely-generated Shift-Invariant (FSI). 

Now, assume that we are given a ladder S :={S^ ■= S^{W 2 ))h>o of SI spaces. Let k > s. We 
say that S provides approximation order k (in 11^2(11'^)), if, for every / G (IR'^), 


dist,(/,5'*):= inf \\f-g\\wURfi ^ 




with the constant C independent of / and h. As is essentially known |16j (and developed fully in 
this paper), the above notion of approximation order depends strongly on k but only mildly on s. 
The ladder S is PSI or FSI if each of its components 5^ is a PSI, or, respectively, FSI space. 

The literature on approximation orders of SI spaces is vast, and it is not within the scope of 
this paper to provide a comprehensive review of it. We refer to the introduction and the references 
of [H] as well as to the exposition and the references in the survey article m Many specific results 
on the topic are reviewed in the body of this article. In particular, a complete characterization 
of the L 2 -approximation orders (i.e., the case s = 0) of PSI ladders is obtained in [Hj, while the 
analogous results for FSI ladders are obtained in P], in and |5]. There are also numerous results 
on approximations in other norms, for example, in Lp. Results and references in this direction can 
be found in [2ni and M- In addition, we refer the reader to m, m, m, isn for information on 
approximation properties of refinable SI spaces, and to HZl, EH], Hi, d, EH, EH] for results on 
wavelet constructions based on SI spaces. 


1.2 Motivation 

While the current level of mathematical understanding of the issue of approximation orders of 
SI spaces is quite advanced, there are numerous gaps and inconsistencies in it. This is exactly 
the motivation behind the present endeavor: obtaining seamless, cohesive (and, so we hope, final) 
theory. We provide a few examples for the “gaps” and “inconsistencies” in the state-of-the-art 
theory. Let us first define two important classes of SI ladders: stationary ladders, and local ones. 

Definition. Let S be an SI ladder. We say that S is stationary if, for every h > 0, = 

S^{-/h) :={f{-/h) : f G S^}. Given a stationary ladder, we say that S is also local if is FSI 
and is generated by a compactly supported 

Discussion 1 (1) Let us assume that S is PSI, stationary and local. Then the entire ladder is 
determined by the (compactly supported) generator of 5^ (since the other spaces in the ladder 
are dilations of S^). In this case, one usually refers to </> as the generator of the ladder. The 
current theory covers the case s > 0, and shows that the approximation order in L 2 (as well as in 
VLI, s > 0) provided by such ladders is intimately related to the order of the zeros cj), the Fourier 
transform of (j), has at the punctured lattice 2tt7Zi'^\0 (cf. §3.2). The smoothness of on the other 
hand, does not play any role, provided, of course, that (/> G L 2 (which is required for the definition 
of L 2 -orders to make sense). Thus, if we replace 4> € L 2 hy its convolution product with a smooth 
generic mollifier, the L 2 -approximation order of the ladder, in general, will not change. In contrast, 
if (j) ^ L 2 while its Fourier transform does have the requisite zeros on 2tt7L'^\0, the smoothing may 
simply result in an L 2 -function, and the ladder may then provide high approximation order in 
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L 2 (despite the fact that the L 2 -approximation order provided by the initial ladder is zero). One 
expects that the extension of the notion of approximation order to W 2 , s < 0 will remove the above 
artificial hump, and this is, indeed, the case. 

(2) Retaining the same setup as in (1), it is also quite well-known that if (pi and (^2 are two 
compactly supported L 2 (lR)-functions, and if the PSI stationary ladder generated by (pj, j = 1 , 2 , 
provides approximation order kj > 0, then the PSI stationary ladder generated by (pi * (p 2 provides 
(at a minimum) approximation order ki + k 2 - One expects then that, if ^2 = 0, the approximation 
order provided by (pi * (p 2 will be at least ki. This, however, is not the case, and there are examples 
when the aforementioned approximation order is smaller than ki. This nuisance is fixed (in §3.4) 
via the introduction of negative approximation orders. 

(3) Let us consider now the case of local stationary FSI ladders in L 2 - In this case = 
S^{L 2 ), with ^ :={(pi,... ,(pr} C L 2 (III‘^) compactly supported, and with = S^{-/h), 0 < h < 
1. A cornerstone in the analysis of the approximation order of such ladders in the existence of 
a superfunction, i.e., a compactly supported function ip & whose associated local stationary 
PSI ladder already provides the same approximation order as the original FSI ladder (cf. §4.2). 
The existence of such a superfunction is proved in [^. However, the Fourier transform of the 
superfunction ip may vanish at the origin, a property that denies us the existence of effective 
numerical approximation schemes from its associated ladder (we refer to such superfunctions as 
“bad”). In [Hj, this problem is overcome, but at the price of imposing an additional condition on 
the vector (its Gramian should be invertible at the origin; cf. H4.5l for a complete discussion). At 
the outset of the current venture, we observed that the condition assumed in [ 3 ] is not necessary for 
the existence of a “good” superfunction (i.e., a superfunction ip for which ip{0) / 0). Unfortunately, 
a good superfunction may not always exist: in ^4.8l we construct an FSI vector (with d = r = 2) for 
which all the compactly supported superfunctions are bad, dashing thereby our hope that a good 
superfunction may be proved to exist in general. 

(4) Let 5 be a ladder as in (3), but assume, in addition, that ;=5'|, is refinable, i.e., that 
S'^ C S^. It is then known (see, e.g., [HI, [13 and m for the PSI case and m for the FSI case) 
that the L 2 -approximation orders provided by the ladder are bounded below by the smoothness 
of if <1> C 11 ^ 2 ( 111 ^)) then the ladder provides approximation order k + 1 or higher. Moreover, 
m proves (for d = 1, and under some mild conditions on < 1 > for d > 1) that approximation order 
/c -|- 1 is implied by the mere existence of a nonzero function f in H ^^(IR'^). However, all these 
results assume more than the smoothness of / and the refinability of 5^: they require in addition 
the entire vector to lie in L 2 - The removal of this condition (§4.9) leads to a conclusion that says, 
essentially, that for S to provide some approximation order, it should contain one nonzero function 
of corresponding smoothness, and nothing else. 

(5) Our final example still deals with refinable ladders. One way to obtain a refinable space 

is to select an r x r matrix P whose entries are trigonometric polynomials and to seek a compactly 
supported vector-valued function <1> that solves the refinement equation <h(2-) = P^. A major 
goal in this direction is to reveal the approximation order of the stationary ladder generated by 
in terms of properties of P (see [HI, [I 3 , 0). The ultimate known result in this direction, 0, 
requires a regularity condition on that necessarily fails once the above refinement equation has 
(in a nontrivial way) more than one solution. Thus, there is no theory at present that deals with 
the approximation orders of refinable vectors, once the refinement equation has multiple solutions. 
§5 deals with the approximation order of stationary refinable ladders and provides a novel theory 
for the case when multiple solutions to the same refinement equation exist. 
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1.3 Layout of this article 

In the introductory m we define the notions of shift-invariance and approximation order and make 
several basic observations that will be extensively used in the sequel. 

Section 121 is devoted to PSI ladders. The section begins with a summary of the known charac¬ 
terization of the L 2 -approximation orders provided by PSI stationary and nonstationary ladders. 
These results are then extended to general W 2 (the end of H3.ll) and connected with the notion 
of the Strang-Fix conditions (EU) and polynomial reproduction f H3.5l) . The results from are 
in turn used in to analyze the dependence of the approximation order notion on the value of 
s, i.e., on the space where the error is measured. The issue of negative approximation orders is 
discussed in mi 

FSI ladders are considered in H It begins, analogously to m with a summary on the L 2 - 
approximation orders of FSI spaces and with the extension of these results to the setting of Sobolev 
spaces. This takes up SH Section focuses on the notion of a superfunction, which is in¬ 
strumental in the reduction of the FSI case to the PSI case. This notion is further used in sa 
to understand polynomial reproduction from FSI spaces and in sa to establish the consistency 
of the notion of approximation order as we operate in different Sobolev spaces. However, not all 
super functions are equally useful, as is made clear in ^4.51 and sa Regardless of whether ‘good’ 
super functions exist in the underlying FSI space, there is an alternative method proposed in sa 
that can always be used to estimate approximation orders. The usefulness of that alternative ap¬ 
proach is demonstrated by an example in sa Section 01 is devoted to refinable FSI spaces. It 
shows that the approximation order of stationary refinable FSI spaces is bounded below by (a weak 
variant of) the decay rate of the Fourier transform of any (nonzero) function in the space. 

In ^ applications of the theory from the preceding sections to multiple solutions to a refinement 
equation are developed. We start by discussing, in HI the structure of the solutions space to 
a refinement question. In HI we introduce the notion of coherent approximation orders, which 
bundles together different solutions to the same refinement equation. In ssini the notion of coherent 
approximation order is associated with a corresponding (novel) notion of universal supervectors; 
those lead to a uniform way of constructing superfunctions in all the FSI spaces that are generated 
by the various solutions to the given refinement equation. 

For the convenience of the reader, some of the notations used in the paper are collected imme¬ 
diately before the bibliography. 

2 SI ladders: the prelude 

We start our analysis with a few elementary, yet very useful, observations concerning the interplay 
between approximation orders in VFI(IR'^) on the one hand, and in L 2 ( 1 R‘^) on the other hand. 

As mentioned before, the symbol IF 2 (IR'^) denotes the Sobolev space of smoothness k. Note 
also the isometry 

: L2(]R'') ^ IF|(IR") : / ^ ((1 + I • . (1) 

Recall that the Sobolev spaces are ordered by embedding: Wl > IT 2 whenever s >t. 
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2.1 Shift-invariance 

The notion of shift-invariance is valid in any fnnction space F whose elements are defined on IR'^, 
and is certainly not specific to ITI Given such a space F, we consider now SI spaces that are 

invariant under integer translations; thus, we refer to a closed subspace S C F as shift-invariant if 
S is invariant under multi-integer shifts 

s e S s(- - a) € S, a € TL'^. 

In agreement with the definitions of PSI and FSI ladders, a principal shift-invariant (PSI) 
space is generated by a single function cp £ F as the closure of 

span[(?!)(- — a) ■. a £ 

in the topology of F, while a finitely generated shift-invariant (FSI) space 5$ is the closure 
of with !> a finite subset of F. 

It is known that an FSI subspace of L 2 (IR‘^) can be characterized on the Fourier domain as 
follows; 

Result 2 (H) For <1 C L 2 (IR.'^), 

S'$(L 2 ( 1 R'^)) = {f £ L 2 (IR'^) : / = r measurable, t{- -|- a) = t, all a £ 2 'k7L'^}. (2) 

That is, the Fourier transform of an element of S'$(L 2 ) is the inner product of two vector-valued 
functions: the vector r (whose entries are measurable and 27r-periodic but otherwise arbitrary), 
and the vector Note that we tacitly assume that the entries of r are indexed by <1 (or by the 
same index set that is used to index <h). 

Since the operators Jg commute with translations, one easily checks that 

= J.sSjML2), (3) 

which, together with Result |2 leads to the following 

Corollary 3 For C IFI 

^^(ITl) = {f £ IT|(IR'^) : / = T*<h, r measurable, r(- + a) = r, all a £ 2 tt7L'^}. (4) 

2.2 Approximation orders 

The basic idea leading to the notion of approximation order is very simple. It is the heuristic un¬ 
derstanding that increasing the density of translations used to define an SI space may improve their 
approximation “power”. At the same time, for numerical reasons (and also for deeper theoretical 
reasons), one would, almost always, change the generator(s) of the SI space when switching from 
to h < I: the new generators should be more localized, and one way, sometime adequate 
sometime not, to modify the generators is by dilation (see the definition of a stationary ladder in 
§ 1 . 2 ). 

The following, simple but important, result connects the approximation orders of SI ladders in 
L 2 to the analogous approximation orders in Wf'- 
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Proposition 4 The ladder S = {S^ = S^{W 2 ))h provides approximation order k in W 2 if and 
only if {JsS^)h provides approximation order k—s in L 2 , where Js is defined as in J7)). 

Proof. Js is an isometry from W 2 to S'® well as from Wf to L 2 . Thus, if JgS provides 

approximation k — s in L 2 then, for every / G W 2 , 

distsif,S’^) = distoiJsf,JsS^)L, < Ch’^-^WJsfW^k-s = ch’^-^Wfiu,. 

'2 2 

Hence 5 provides approximation order k in Wf. The converse is proved in the same manner. □ 

As already indicated before, the two most important cases of SI ladders are 

• PSI: each is an /i-dilate of some PSI space, i.e., = S^^{-/h); a PSI ladder may be 

stationary or nonstationary depending on whether or not the generator of is independent 
of h. 

• FSI: each is an /i-dilate of some FSI space 5$^; an FSI ladder, just like a PSI ladder, may 
be stationary or nonstationary. 

Nonstationary FSI ladders are broad enough to cover almost all situations of interest in ap¬ 
plications. Thus it is of primary importance to be able to characterize the approximation orders 
provided by such ladders. It turns out that nonstationary ladders are useful not only on their own, 
but also as a tool for analyzing stationary ladders. 

Corollary 5 An FSI ladder {S^ := S^^{-/h))h provides approximation order k in Wf, if and only 
if the FSI ladder {S^f^{-/h))h, T/j :=(1 -|- | ■ provides approximation order k—s in L 2 . 

Proof In view of Propositionwe only need to identify {JsS^){h-) as with T/j as in the 
corollary. Now, by Corollary 01 / G iff / € ^2 r 27r-periodic. Thus, 

/ G 5^ iff / G Wf and / = T*^h{h-), with r 27r//i-periodic. Thus / G JsS^ iff / £ L 2 , and 

f = {l + \-\y/\*$h{h-). 

Dilating the last equation, we obtain that / G {JsS’^){h-) iff / G L 2 and 

f = {i + \-/hfy/\*^h 

for a 27r-periodic r. By Result^ this last condition is equivalent to / being in 5ii(^(L2)- D 

Note that the ladder associated with {^h)h in the above result is nonstationary even when we 
assume the original one to be stationary, i.e., when we assume to be independent of h. 

3 PSI ladders 

We start our study of PSI ladders by recalling the characterization of the L 2 -approximation order 
of these spaces. We then extend the result to the Sobolev space Wf- The general result is then 
connected with the notions of the Strang-Fix conditions and polynomial reproduction. In turn, 
those latter notions allow us to understand the dependence of the approximation order notion on 
the value of s, i.e., on the space where the error is measured. 
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3.1 Approximation orders of PSI ladders 

Note that the first part of the next result is not entirely a special case of second part (although it 
can be derived from it with ease). 


Result 6 m Theorems 1.6 and 4.3]) 

1. The stationary PSI ladder S = {S^ := {L 2 )), with = S^{-/h), cj) & 1^2(51'^); provides 

approximation order k if and only if there exists a neighborhood H of 0 sueh that 


[^, 0 ] 



G Lcaiyi)- 


Here [^>] := l^(-+ a)P, l^(-+ «)P- 


2. The nonstationary PSI ladder S = {S^ := S^{L 2 )), with S^ = S^^^{-/h), 4>h £ 1^2(111'^); P'^o- 
vides approximation order k if and only if, for some /iq > 0 and some neighborhood 11 of 0, 
the collection of functions 


1 


0 < h < ho, 


lies in Loo(ll) o,nd is hounded there. 


Combining Proposition 0 and Result El we obtain the analogous result for Sobolev spaces. 


Theorem 7 Let s G IR and k > s. Assume also that k is nonnegative. 

1. The stationary PSI ladder S = (5^ := ))? — S,j,{-/h), cj) G 11^2^, provides ap¬ 

proximation order k if and only if there exists a neighborhood H of 0 such that 

:= ^321 e Loo(l!). (5) 

Here [^, $]s := Eag 27 rffi'* l^(' + «)Pl ' +«P^ E„g2,rffi‘*\0 l^(- + “)l^l 

2. The nonstationary PSI ladder S = {S^ := S^{Wf)), with S^ = S^,^{-/h), 4>h £ IPI; provides 
approximation order k if and only if for some h^ > D and some neighborhood H of 0, the 
colleetion of funetions 


[4>h, _ 1 

^hAh]s,h (I • P + 


0 < h < Lq 


lies in Loo(ll) and is hounded there. Here, [(fh, <kh]s,h '■=[(t>h, (fh]^s + I'/’/iPd ‘ P + 


( 6 ) 


Proof. The second part of the current theorem follows from the second part of Result El and the 
PSI case of Corollary El Together, these two results yield the requisite characterization, but with 
[4>h,4>h]^ replaced by Z]ag 27 rS‘^\o + «)P(| • +ap + h'^Y. However, for a 7 ^ 0 , we can replace 
(I ■ Tap + h'^)^ by its equivalent expression | • +ap®. 










It remains to show that in the stationary case, i.e., when (jjh = (j) for all h, is equivalent 
to 0. The fact that the former implies the latter is obvious (one simply should take /i ^ 0 in 
0 and invoke the uniform boundedness of the collection of functions that appears there). For the 
converse, we observe that (when (j)h := (p for all h) the uniform boundedness of the functions in 0 
is equivalent to the validity of the inequalities 

|^|2 - + 

for some positive c > 0. Moreover, since we assume k — s > 0, we can force (1-1^ + < c by 

making h small enough and changing if necessary. This leaves us with 

< C(| • |2 + /i2)^ a.e. 

w ~ 

as the requisite boundedness. This is definitely implied by ©, as the left hand side in the display 
above is independent of h and since A: > 0. □ 

Remark on notation. For brevity, we will use in the sequel the expressions provides 

approximation order k' and ‘0 provides approximation order k in ITI’ to mean that the stationary 
ladder generated by provides approximation order k in W^- 

3.2 Strang-Fix conditions 

Given 4> G VFKlR'^), and k > 0, one says that p satisfies the Strang-Fix (SF) condition 
of order k f|36ji. if p has a zero of order k at each point a G 27rK'^\0. It is well known that 
the L 2 -approximation order of a stationary PSI ladder is closely related to the order of the SF 
condition satisfied by the generator (p of the ladder. To be precise, a full characterization requires 
a nondegeneracy condition on p at the origin. First, let us cite the L 2 -result. 

Result 8(0 Theorems 1.14, 5.14]) Assume that 0 < rj < \<p\ a.e. on some neighborhood 
of the origin. Let d := + 27rffi‘^\0. If (p € IF 2 (^) for some p > k + d/2, then S^{L 2 ) provides 

approximation order k (in L 2 ) if and only if (p satisfies the Strang-Fix conditions of order k, i.e., 
near the origin 

\P{. + a)\ = 0{\-f) for all a G 27rS'^\0. (7) 

Here W^iA) is the local version of IT2(1R'^); cf. P Chapter 7]. For our purposes, it is only 
important that the norm on W^iA) has a subadditivity property, i.e., 

X! \\f\\wf(a+n) - const (8) 

a 

and that the Sobolev embedding theorem for such spaces still holds, in particular, that the bounded 
(compact) embedding 

W^{a + n)^W^{a + n) (9) 

is valid. Note that the condition p G W^iA) is weaker than the more traditional decay condition 
on p 

\P\ = 0{\-\-^-'^-^), e>0. 
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which implies global smoothness of cj). 

We now show that the Strang-Fix conditions also characterize approximation power in a Sobolev 
space. 


Theorem 9 Let k > 0, s < k, (j) € W^- Suppose that, for some 771,772 > 0 and for some ball Ll 
centered at the origin, 

7?i < \4>\ < ^?2 a.e. on Q (10) 

U\\l,A-= (11) 

where A denotes the set Ll + 27rffi'^\0 and denotes the monomial derivative of order 7 G 
normalized by 7 !. Then provides approximation order k (in Wf) if and only if 0 holds. 


Proof. Set 

R:=\-\-^^ ^(■ + /3)P|-+/3p. (12) 

/7e27rS‘*\0 


Suppose (p provides approximation order k in Wf. Then © 
a.e. on fl, 


R 

W + R 


0(1 • |2^-2"). 


holds by Theorem [3 or equivalently, 


Since k > s, and (p is bounded on Q, we conclude that, around the origin, R = 0(| • ^^). This 

readily implies ©• 

Now suppose (p satisfies ©. With R as above, we invoke m to conclude that 


m\ 


Loo 


<o|- 


|2A:—2s I 


\k,A 


= 0(1 


\2k—2s\ 


However, the left-hand side of © equals 


I , ^2s—2kj^ 

|0|2 + i2 ■ 

We have just argued that the numerator in this expression is bounded. The denominator of the 
expression is bounded away of zero thanks to m- This implies ®. □ 

Note that condition m was required only for the “if” implication in the above result. 


Corollary 10 In the notation of Theorem\^ let p > k + d/2 and let H := H -|- 27rK'^\0. Then the 
conclusions of Theorem\^ remain valid when we replace condition m by: 

(i) for s < 0, the condition that (p G IT 2 (A). 

(a) for s >0, the condition that (1 -|- | • p)^/^0 G ITI(IR*^), or the stronger condition that <p G Wf 
and (p = 0(1 • 1 “*^“'^“^), e > 0 , at 00 . 
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Note that the first condition in (ii) above implies, whenever s > 0, that (j) € hence is 

stronger than the condition assumed in (i). 

Proof. It is clearly sufficient to prove that each of the conditions in (i) and (ii) implies CU). 

(i) : Using Q, together with the fact that that the sets (/3 + fl) are all translates of 0, the 
Sobolev embedding theorem applies to yield that, for s < 0, 

\\$\\k,A^Ci X! \(^\‘^'^\\$\\wPin+i3) - X! \\$\\wPin+i3)- 

/3G27rffi‘*\0 /3G27r®‘*\0 

The right-hand-side in the above is bounded, thanks to (jHl), by a constant multiple of ||</’||^/=(^p 
Hence condition m is satisfied. 

(ii) ; The second condition in (ii) clearly implies the first one. Now assume the hrst condition 

in (ii), i.e., that / :=(1 -|- | • G lUI• Then / is locally in L 2 and 

l/^P^II/lli2(/3-i-20) ^ C'll/llw| < °°' 

/3e27rffi'^\0 

However, ||<('||u/'’(/3+n) — C'||/||l 2 (/ 3 + 20 )) the argument in the proof of (i) then applies to yield 

(fTT|). " □ 

3.3 Approximation orders are independent of the underlying VU 2 ^ space 

We are now in a position to observe that the definitions of approximation order, if made with 
respect to different Sobolev spaces, are consistent in the following sense. 

Proposition 11 If S^iW^) provides approximation order k > 0, k > s, then S^{W 2 ) provides the 
same approximation order for any t < s. 


Proof. First note that (f is an element of whenever t < s, since Wf is embedded in W 2 . 
Now note that, by Theorem [7| (f provides approximation order k if and only if ((21) holds. The 
left-hand side of © satisfies 

(1 - M^,sl • E l^(- + • P"- (13) 

/3e27rffi‘^\0 

Since ^ const,^^^, the set fl can be assumed to be small enough so that, e.g., 

1 — M.cj,^s\ ■ > 1/2 a.e. on H. 


Then 

E \k- + /3)Pl • +/3|'* < E l^(- + ^ 

/3e27rffi‘*\0 /9G27r®‘^\0 


This implies 


< 


S/3g27r®‘^\0 




\2k 


< 2M^ 


,s ■ 


Thus, 5*0 provides approximation order k also in W 2 . 


□ 
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Proposition im shows that (p provides approximation order on the whole interval {VPl ■ t ^ s} 
of Sobolev spaces once it does so in the space W^- 

Let us now show that, under the regularity assumptions already used in Theorem [7J a converse 
also holds. 

Theorem 12 Let t < s < k, k > 0. Suppose cp G VFI and that it satisfies (with respect 

to s). Then SfjyiW^) provides approximation order k iff SfjyiWfi) provides that same approximation 
order. 

Proof. The “only if” implication was proved in Proposition II II without appealing to (ll()l) - (llll) . 
We prove the “if” assertion as follows. First, since p provides approximation order k in VFl, 
while satisfying m, it must satisfy the SF conditions or order k (we do not need m for that 
part), by virtue of TheoremlHl Then, once p satisfies the SF conditions of order k, the facts that it 
belongs to Wl and satisfies (fTU)) -(ITT]) imply, again by Theorem that it provides approximation 
order k in Wf. □ 

3.4 Negative approximation orders 

When would it make sense to have a shift-invariant space that provides negative approximation 
order? 

Suppose we form a convolution of two compactly supported distributions pi, pi{0) 0, i = 1,2. 

If pi, i = 1,2, provides a positive approximation order ki > 0, then their convolution p provides 
approximation order ki + k 2 . This is to be expected if one assumes that the Strang-Fix conditions 
are equivalent to approximation power (which is almost true): then p = pip 2 satisfies the Strang- 
Fix conditions of order ki + k 2 whenever each pi satisfies the Strang-Fix conditions of order ki. A 
sample of a rigorous statement in this direction is as follows: 

Proposition 13 Let pi € Wfp i = 1,2 be nondegenerate compactly supported distributions, i.e., 
pi{P) p 0, and provide approximation order ki > 0, ki > Si, i = 1,2 in their respective spaces. 
Then p:=pi * p 2 provides approximation order ki + k 2 in . 

Proof. First observe that p G Indeed, 

< iii?ii(i +1 ■ i")*‘-'"iuiii&i{i +1 ■ (n) 

Since the pfs are compactly supported, their Fourier transforms are entire functions. Moreover, 
the products pfil + \ ■ |2)^i/2^ i = 1,2, are in L 2 , and, therefore, their inverse transforms are in 
L 2 , too. Those inverse transforms are the result of applying a singular convolution operator to p. 
Since the convolutor decays rapidly at 00 , and since p is compactly supported, the result decays 
rapidly at 00 . Altogether, we conclude that {pi{l -|- | • is in Li, and consequently each 

pi{l + I • must tend to zero at infinity. Therefore their Loo-norms must be finite. So, the 

right-hand side of m is finite, hence p is in . 

Now, since we assume pi to provide approximation order ki, and since pi is bounded around 
the origin, then (cf. the first part in the proof of Theorem [HI) 

= i = l,2. 
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(15) 


where we recall the notation [ 5 , 5 ]^ '■={g-,9\s — IffPl • P* from Theorem |71 Bnt 


hence 

Invoking the fact that '0(0) 7 ^ 0, we finally conclude that 

[0; '^]si+S2 _ 0(\ ■ 2(si+S2)'\ 

[ip,^]si+s2 

This, in view of Theorem [7| finishes the proof. □ 


Now, what if, upon convolving a given distribution 0i with another distribution 4>2 one discovers 
that the approximation order of 0i * 02 is smaller than that of 0i? Then it is natural to assign a 
negative approximation order to the distribution 02 - It makes little sense to define the notion of 
negative approximation order in terms of the ability to approximate functions. We choose, instead, 
the following technical definition, which is consistent with the discussion so far, as well as with the 
argument used in the proof of our last result. 

Definition. Let s < k < 0, and let 0 G W^- We say that the (stationary ladder generated by) 
0 provides approximation order k (in W 2 ), if, for some neighborhood of the origin. 


M 


1041 

[05 0]s 


1 

2k—2s 


G Loo(n). 


Note that the above definition is consistent with the case A: > 0. In this case, the definition of 
approximation order is different, but the characterization provided in © of Theorem [7| is exactly 
in the same terms. 

Equipped with this last definition, we can extend Proposition^] as follows; 

Proposition 14 Let (jii G IT^^IR'^), i = 1,2 provide approximation order ki > Si in IT^^IR'^), 
i = 1,2. If the convolution product 0 := 0i * 02 lies in 1 T 2 ^"*~^^(]R'^), then it provides approximation 
order ki + k 2 there. 


Proof. Since ki > Si for i 
order that 


1 , 2 , it follows directly from the extended definition of approximation 


[0i5 0i]° ^ Qn |2fci 


i = 1,2. 


These two estimates, together with the inequality m, imply that 




(90 . | 2 fci+ 2 fc 2 ^^ hence 


[ 0 > 0 ]° 1 +S 2 


| 0 | 2 | . |2si+2s2 _j_ [ 0 ^ 0 ]O 


< 


S1+S2 


(/;|2] . |2si+2S2 


= o(|. 


|2/ciH-2Ai 2—2si—2 s2 > 


This completes the proof. 


□ 
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Corollary 15 Let cjii G i = 1,2 he compactly supported distributions and provide ap¬ 

proximation order ki > Si in iy2*(lR'^), i = 1,2. Then (/>i * 4>2 provides approximation order ki + /c2 
in 

Proof. This fact follows from Proposition 1141 since we know already from the proof of Proposi¬ 
tion cni that the convolution of two compactly supported distributions in ^ = Ij 2 , lies in 

Remark. In the rest of the paper, we only consider, by default, generators </> of stationary 
ladders that provide approximation order no smaller than 0. Note that this is the case when (p is 
of compact support and satisfies 4>{0) 7 ^ 0 . 

3.5 Polynomial reproduction 

We restrict our attention in this subsection to local stationary PSI ladders, and focus on the 
properties of the compactly supported generator cp of the underlying ladder. To be sure, all the 
results here extend, almost verbatim, to generators (p with sufficient decay at 00 , for example, 

= 0{\ ■ at 00 , with k the investigated approximation order and e > 0 . 

The theory of approximation orders of local stationary PSI ladders focuses, and rightly so, on 
the satisfaction of the SF conditions (cf. il3.2[ and also the application of those conditions in . 
Under the compact support assumption on p, the SF conditions are known to be equivalent to the 
polynomial reproduction property, the latter being the subject of the current subsection. ^ 

The connection between the SF conditions and polynomial reproduction is classically known, 
and can be dated back to Schoenberg (d = 1, [22] ), and Strang and Fix, [22] • See also |21. Our 
approach here follows [7j. Altogether, the results of this subsection are included for completeness, 
especially since the polynomial reproduction property in the PSI case is key to the understanding 
of the more complicated polynomial reproduction property of FSI spaces (S31), as well as the sum 
rules of refinable FSI spaces i JOl . 

Suppose that p is compactly supported. Let us first attempt to connect the approximation 
orders provided by its stationary PSI ladder to the SF conditions. To this end, we would like 
to invoke Theorem El This theorem requires the satisfaction of m and m- Condition m 
is satisfied once p G W2, as (ii) of Corollary ^1 shows. The fact that a compactly supported 
distribution belongs to some W2 is well-known, and follows from the fact that it is necessarily 
of finite order (as a distribution). As to (llOj) . since p is continuous, this condition is presently 
equivalent to the nondegeneracy requirement 

^( 0 ) + 0 . 


Thus we obtain the following result: 

Corollary 16 Let p be a compactly supported distribution, and assume that p{0) 7 ^ 0. Then there 
exists s G IR such that p G Wl. Moreover, the following conditions are then equivalent, for any 
given k > 0: 

^Prior to the publication of and [3|, approximation orders of stationary PSI ladders were usually derived directly 
from the polynomial reproduction property, while the SF conditions were considered to be a technical way for the 
verification of polynomial reproduction. However, as the discussion in this article clearly shows, the SF conditions 
characterize the approximation orders of the ladder even when a slow decay of the generator p renders the polynomial 
reproduction property meaningless. 
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(i) (j) satisfies the SF conditions of order k. 

(ii) The stationary PSI ladder generated by cj) provides approximation order k (in Wf)- 
Now, recall that reproducing polynomials of total degree less than k means that 

4> *' n<fc c n<fc. 

The symbol *' denotes the semi-discrete convolution 

9*' ■ f ^ 9{- (16) 

n := n(]R'^) is the space of all d-variate polynomials, and n<fc :={p G II : degp < k}. 

One way to connect the polynomial reproduction to the SF condition is via the following variant 
of Poisson’s summation formula, m, 

4>*'f= e„ : X (17) 

which is valid for every compactly supported f and every C'°°-function / (the convergence of the 
right-hand-side series is in the topology of tempered distributions). Now, for a polynomial /, one 
easily verifies that (j) * (e-af) = 0 iff has a zero of order deg / -|- 1 at a. Thus, once (j) satisfies 
the SF conditions of order k, we have that (f *' f = (f) * f for all / G n<fc. This establishes the 
sufficiency of the SF conditions, since (p* always maps into itself. On the other hand, ii cp*' f 
is a polynomial of degree < fc, then (dzi) shows that 

4>*ieaf) (18) 

is also a polynomial of degree < k. This is possible [3 Proof of (2.10) Lemma] only if all the 
summands in m vanish. In conclusion, cp satisfies the SF conditions of order k if and only if 

(pfi = (p*^ on n<fc. 

Since, as we already said, is an invariant subspace of (p* (with or without the SF conditions), 
we finally need only to guarantee that (p* be injective on polynomials, or equivalently, we need to 
assume that ())(0) 0. Indeed, the condition (p{0) 7 ^ 0 is necessary and sufficient for cp* to be an 

automorphism on n<fc (for any positive integer k), and we arrived at: 

Theorem 17 Let (p be any compactly supported distribution with (p{0) 0, and let k be a positive 

integer. Then cp provides approximation order k in some Sobolev space W 2 , s <k, if and only if it 
reproduces polynomials of total degree less than k. 

Remark. As alluded to before in a footnote. Theorem ITTl could also be proved directly, avoiding 
the use of the Strang-Fix conditions and constructing instead a quasiinterpolant Q : Wf S{(p) 
such that Sp = p for any p G for a detailed discussion of this method see [3 Section 4]. 
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4 FSI ladders 


We start this section, just like in the PSI case, by recalling the characterization of the L 2 - 
approximation orders of FSI spaces and extending the result to the setting of Sobolev spaces. 
We then focus in 114.2l ou the notion of a superfunction, which leads to the reduction of the FSI case 
to the PSI case. Besides, this notion proves to be very helpful in understanding polynomial repro¬ 
duction from FSI spaces (see EH). In our setting of W 2 , it also helps to establish the consistency 
of the notion of approximation order as we vary s (see ESI). However, not every super function can 
be used for these and/or for other purposes, and this brings one to the notions of ‘good’ and ‘bad’ 
superfunctions that are discussed in ESI We show, in ESI that there exist FSI spaces that do not 
contain any good superfunctions. Regardless of whether or not good superfunctions are around, 
an alternative approach, that is presented in ESI can always be used to bound the approxima¬ 
tion order from above. The efficacy of this method is demonstrated in ESI where we recover the 
well-known example of C^-cubics on a three-directional mesh, [0]: this is an example of a bivariate 
stationary local FSI ladder that, while reproducing all polynomials in n 3 , fails to provide the “ex¬ 
pected” approximation order 4. Finally, H4. 91 applies the results obtained in this section to the case 
when the vector is refinable: in establishes a lower bound on the approximation order provided 
by in terms of the decay of the Fourier transform of any non-zero function in S^. 

4.1 Characterization of approximation power 

The first three results of this section form a summary of the known characterization of approxima¬ 
tion power valid in L 2 , while the rest constitutes the characterization in the more general setting 
of IF|. 

Result 18 (|Sl Theorem 2.2]) The stationary FSI ladder S = {S^ := S^{L 2 )), with = S^{-/h), 
$ C L 2 (IR‘^); provides approximation order k if and only if there exists a neighborhood H o/O such 
that 

GLoo(f^). 

Here 

G$:= ^ + + = (19) 

Also, the expression is taken to mean any solution to the equation = $. 4 simple 

linear-algebraic argument shows that the latter equation is always solvable whether or not G# is 
invertible, since one of the rank-one terms in GB is <^4**. 

Result 19 (O Theorem 2.7]) An FSI nonstationary ladder S = {S^ := S^{L 2 )), with S^ = 
S^hi'/h), C L 2 (IR‘^)x provides approximation order k if and only if, for some /io > 0 and some 
neighborhood H of the origin, the collection of functions 

(1 - (I . |2 _|_ ’ h < ho 

lies in Loo(H) and is bounded there. 
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We also require the following equivalent formulation of the last characterization, in which we 
use the notation: 

(70:= 5(-+ «)$*(• +a) ( 20 ) 

oG27rffi'^\0 


Result 20 (another version of B£sult \19\) The FSI nonstationary ladder S = {S^ := {L 2 )), with 

= S^^{-/h), C L 2 (IR‘^); provides approximation order k if and only if the collection of 
functions ■ 0 < h < ho), where 


■M.^,s,h '■ I—> 


1 . p v*G%{u)v 

(Iwp + ijec* v*G^{u})v 


is bounded in -Loo(^ 2 ) for some neighborhood Q of the origin and some ho > 0 . 


Using these results, one obtains the following characterization of approximation power in Wf. 

Theorem 21 

1. An FSI stationary ladder S = :=S’^{W 2 )), with = Sq,{-/h), C Wf, provides approx¬ 

imation order k > 0 if and only if there exists a neighborhood 0 0 / 0 such that the function 

1 J(w)v 

: a; inf \ , belongs to Loo(f 2 ). ( 21 ) 

zs v*G^^s[^)v 

Here 

G*<i>,s ■= $(•+ «)$*(• +a)| •+ap^, 

QE27rS'^ 

C'L •= E $(- + a)$*(- + «)l 

ae27r®‘*\0 

2. An FSI nonstationary ladder S = {S^ := S^{W 2 )), with = S<i,^{-/h), C Wl, provides 

approximation order k > 0 if and only if there exists a neighborhood H of 0 and ho > 0 such 
that the collection of functions : 0 < h < ho), with 

1 v*G% u{u;)v 

: tv ^ ^ inf , is bounded in Loo(f 2 ). ( 22 ) 


Here, 


G$h,s,h(‘^) 


G'$^,s,h(<^) 


^>/i(w + a)$/i*(a; + a)(|a; + ap + h^)^ 
Y ‘h/j(a; + Q;)<h/j (a; + a)(|a; + ap +/i^)"". 

ae27ra‘*\0 
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Proof. The proof is analogous to that of Theorem[71 In particular, part 2 of the current theorem 
is a direct consequence of Result EDI and Propositional 

Now we use the result of part 2 to derive part 1. In the stationary case, = <1> for all h, so 
the left-hand side of (|^ becomes 

+ ^ n* + “)(l^ + “P + ^ 

Since the numerator of the infimum expression is bounded above and below by positive multiples 
of gV, the collection is bounded in Loo(^^) if and only if the collection of functions 


U) h-!■ 


1 _ 

(|wp -|- V v*G^ ^(^uj)v + + h‘^)^v 


(23) 


is bounded in Loo(f^)- Since k and k — s are nonnegative, for a fixed v and a fixed u G fl (assuming 
fl is sufficiently small), the expression 


1 _ ^*G|,,(c^)^ _ 

(Iwp + /i2)fc-s u*G|_^(n;)n + n*($$*)(a;)(|a;|2 K^Yv 


monotonically increases (as /i —> 0) to 


1 n*G^_^(a;)n 

|^| 2 (fc-s) y*G^^s{^)v' 

hence (Eni monotonically increases to the function in (ED). Therefore, the collection (A^</,,s,fe) is 
bounded if and only if m holds. □ 

Remark on notation. As in the PSI case, we shall use in the sequel language such as ‘an FSI 
space S'$(IT 2 *) provides approximation order k' and even ‘<1> C IT 2 * provides approximation order k' 
and will mean by that that the FSI stationary ladder generated by ^^(Wl) provides approximation 
order k in W^- 

The L 2 -characterizations above (Results [THIEn|) are connected to the notion of superfunctions. 
We will now discuss this notion, and extend it to the setting of Sobolev spaces. 

4.2 Superfunctions 

Let S be an SI space, and let S = {S^ :=S(-/h)) be the associated stationary ladder. A function 
g G S C L 2 is a superfunction in S if the PSI stationary ladder it generates provides the same 
approximation order as that of S (or, more precisely, of S). For sure, L 2 in this definition can be 
replaced by any Sobolev space 

The question of existence of superfunctions in FSI spaces can be answered in the affirmative 
using Theorem 1211 

Theorem 22 Any FSI space 5# C LFKIR®*) contains a superfunction. 
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Proof. Let 0 be as in Theorem inn For each fixed a; G there exists a vector t'o(w) G C*" of 
(e.g., Euclidean) norm 1 that minimizes the ratio v*G^ g{uj)v/v*G^^s{<.^)v- Now extend vq to the 
cube [—TTjVr)'^ in an arbitrary way, provided the norm uq is everywhere equal to 1. Finally, extend 
vq so defined to a 27r-periodic vector-valued function. 

Now, suppose that provides approximation order k. Then, in view of Theorem 121 L the vector 
Vq satisfies 

< const | • -I- ■ P®) a.e. in Q, 

or, equivalently 

(1 — const I • < const | • p^UQ‘f>^>*uo a.e. in fl. 

By changing fl if needed (and using the fact that k > s), we obtain that 

< G\ ■ a.e. in Q 

for some constant G. This implies that, for almost every fixed a; G O, the smallest eigenvalue of 
the measurable Hermitian matrix H{uj) :=(G^ ^ — G\- is nonpositive. By Lemma 2.3.5 

from m, we can define a map rc on 12 such that (i) for almost every a; G 12, w{uj) is a normalized 
eigenvector of H{uj) that corresponds to the minimal eigenvalue, and (ii) w is measurable on 12. 
Without loss, we assume that our original vq coincides with w on 12. In particular, vq is now known 
to be measurable. 

Let (j) be the (scalar) distribution whose Fourier transform satisfies <j) = Vq^. To show that cj) is 
a superfunction for 5$, we only need to verify that it belongs to 5$, since it follows directly from 
the construction of cj) that the space provides approximation order k. Since ^ = Ug<I>, we only 
need, in view of Corollary to show that 4 > G W^- This final result is a simple consequence of the 
representation (/> = using the facts that C W2 and that the entries of vq are bounded. □ 

This theorem extends the known L2 result. However, in the L 2 -case, a superfunction was 
originally constructed as the orthogonal projection : L 2 —> of the sinc-function 


d 

sinc(x) := 

i=l 


sin(7rx(f)) 

Trx{i) 


The fact that Pi>(sinc) is a superfunction follows from the general principle: 


Result 23 m) Let S^{L 2 ) be an FSI space that provides approximation order ki > 0, and let 
Sg{L 2 ) he a PSI space that provides approximation order k 2 > 0- Then the PSI space generated by 
the orthogonal projection P^g on S'# in L 2 provides approximation order m.\n{ki,k 2 }. 


If, in Result above, we choose g such that k 2 > ki, and ^#(^ 2 ) does not provide an ap¬ 
proximation order greater than ki, then P^g is a super function. It is easily checked that the 
approximation order provided by the space Ssinc is infinite (i.e., exceeds any finite /c), hence the 
following corollary. 


Result 24 m) An FSI space S#(L 2 ) provides the same approximation order as the PSI space 
generated by P#(sinc) in L 2 . 

To find the L 2 -projection on a FSI space of a function /, one solves equation (|M . 
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Result 25 ( 0 ) The L 2 -projection P<^{f) of f G L 2 on an FSI space satisfies 




( 24 ) 


with Tf any solution of 

Gq>Tf = [$,/]. 

Here, Tf is a vector-valued function (indexed hy whose entries are measurable and 2 t:- periodic, 
and the symbol [$,/] stands for ([0,/])</,£$, where [f ,g\ ■.= Y.a& 7 L'^ /(' + 2 ™) 5 (- + 27ra). 

Remark: Results 1221HH and 1221 are all corollaries to Theorem 3.3 of [21. 

The above results extend easily to Sobolev spaces. Indeed, Results OHl and 121 require only one 
assumption, viz. 

PaPq = PpAgPg, (25) 

where A denotes an arbitrary SI subspace of Wf, and Pa, Pg, Pp^g are the orthogonal projectors 
from IT 2 ^ onto A, Sg{Wf), Sp^gfWf), respectively. 

Under this condition, the analysis from j21 Section 3] leading to Results (21 and |21 goes through 
verbatim. Since the lUI-version of ((221) is a simple consequence of the L 2 -version when combined 
with the identity Pj^aJi = JiPa (where A an SI subspace of Wf and Pj^a the orthogonal projector 
onto the space JtA in L 2 ), we obtain the following extension. 

Theorem 26 Let ^^(VUI) ^0 an FSI space that provides approximation order ki > 0, and let 
SgfWf) be a PSI space that provides approximation order k 2 > 0. Set fi:=P,^g, with Rj> the 
orthogonal projection ofWf onto S#(IU|). Then the stationary PSI ladder generated by provides 
approximation order m.\\i{ki,k 2 }. Specifically, for every f G Wf and h> 0, 

disU(/,5^(IT|)) < disU(/,5|(IU|)) + 2disU(/,5,^(IU|)). 

In particular, ^^(ITI) provides the same approximation order as the PSI space generated 6 ?/P$(sinc) 
in Wf. 

Superfunctions are obviously useful if one wishes to approximate functions from a given SI space 
S, for if a superfunction is known explicitly, one can instead approximate from the simpler space 
Sf,. In addition, it is already well established in the L 2 -theory that superfunctions give rise to quasi- 
interpolants, i.e., bounded linear maps into the underlying SI space S that reproduce polynomials 
contained in S (see, e.g., 0, 0). Superfunctions were also used in |21, |21 as purely theoretical 
tools. 

The natural expectation is that superfunctions play similar roles in the setting of Sobolev spaces. 
That turns out to be the case. In particular, the superfunction method allows us to lift painlessly 
various results from the PSI setup to the FSI one. This includes the discussion concerning the 
consistency of the definitions of approximation orders in different Sobolev spaces, which we embark 
on in the next subsection. 

4.3 Approximation orders are independent of the underlying Wf space 

Proposition 27 If an FSI space S'$(IT 2 ) provides approximation order k > 0, k > s, then S^fW^) 
does so for any t < s. 
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Proof. Assuming that provides approximation order k, Theorem 1221 ascertains that 

contains a PSI subspace that already provides approximation order k. Since 

P’ easily follows from Corollary |31) , </> G too. By Proposition ITTl 

S^iW^) then provides approximation order k in W 2 , therefore S^iW^) provides approximation 
order (at least) k in W^- □ 

Similarly to the PSI case, a converse also holds under some regularity assumptions on the 
super function. 

Proposition 28 Let k > s > t, k > t). Suppose (j),^ C 11^2^. Suppose that provides 

approximation order k in W 2 , and that </> G is a corresponding superfunction. If (p satis¬ 
fies then 5$(IT2*) provides approximation order k in Wf, too. 

Proof. We apply Theorem 1121 to the function p to show that provides approximation 

order k, which implies that 5$( 1 ^ 2 ^) also provides approximation order (at least) k. □ 

The theorem highlights a central point: it is useful to know that an FSI space contains a “good” 
super function. In the current context “good” in interpreted as “satisfying (ITni)- (|TT]) ”. We will come 
back to this issue later, but hrst we show how the superfunction method reduces the polynomial 
reproduction issue in the FSI setup back to the simpler PSI setup. 

4.4 Polynomial reproduction 

Let be a vector of compactly supported elements in Wf, s G IR. Suppose that <I> provides approx¬ 
imation order A: > 0 in Wf. Let us assume, further, that 5(i>(IF2®) contains a good superfunction ip 
is the sense that: 

1. p^O) p 0, and 

2. '0 is a finite linear combination of the shifts of 4 (hence, in particular, is compactly supported). 

We note that the current notion of “good” is stronger (i.e., implies) the one that was discussed at 
the end of the last subsection (as the argument in 51 showsL 

By our assumptions here, p = v*^, with v a vector of trigonometric polynomials. Therefore, 
with (a(^)(^g$ the Fourier coefficients of the entries of v, we have the representation 


P = P *' a<p, 

and each Oip : —> C is finitely supported. Here, *' is the semi-discrete convolution, (|T?)1) . 

Next, since p{fS) p 0, and *S'^(IF 2 ) provides approximation order k, we conclude from Theo¬ 
rem im that p*' maps H^^ onto itself. Writing p*' in terms of we obtains 

P *' f = ^ *' ^<1^^ *' f ^ ^^<t> *' /)• 

The above representation leads to several conclusions that we summarize in our next result: 
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Corollary 29 Let ^ be a compactly supported vector that provides approximation order k > 0 in 
W 2 , and assume that S'$(VF 2 ) contains a good superfunction in the above sense. Then there exist 
finitely supported sequences a^j, : ^ C, (/> G $ such that, for every f G 

is a polynomial. Here, a^* f denotes the discrete convolution of atj, and fi^,^- 
The polynomial Tf is identical to the result of the following continuous convolution 

J2f>*ia^*f) = J2 * /)(• - j)- ( 26 ) 

Moreover, the map Tk ■=T\ is an automorphism. 

There are several immediate conclusions that can be derived directly from the above corollary. 
For example, since T^. can be extended to a convolution operator, it commutes with differentiation 
in the sense that W^T = TD^ for every 7 G and commutes also with translations. 

A simpler consequence is as follows: since is an automorphism, then every monomial ()“, 
|q;| < k, lies in its range. (Here, the symbol ()“ stands for the normalized monomial 

0" (x)“:=x“/a!. 

We also use in the sequel for the normalized monomial derivative.) Thus, the following is true: 

Corollary 30 Let <1> be as in Corollarv\2fA Then there exist polynomials {ga)a£'Z’‘ such that, for 
|q;| < k, 

^ {a^*ga) = ()“• 

The result shows that every ()“ is writable as fi*' f<j>,a-, for suitable polynomials {f(i,^a)(j>- How¬ 
ever, the result shows more: it decomposes each into ga, with the first factor independent 
of a (and is finitely supported), and the second independent of cf (and is a polynomial). 

The reader might wonder how realistic the assumption about the existence of a good super¬ 
function is. We discuss that issue in this section as well as in §5. A sufficient condition for the 
existence of a good superfunction as above is the invertibility, in a suitable sense, of the Gramian 
around the origin. We also note that our results here recover the results of m Section 3] (cf. 
also i, nni). The underlying assumption in m is that the shifts of the distributions (/> G are 
linearly independent, a condition that is significantly stronger than the Gramian invertibility that 
we have alluded to above. At the same time, our derivation here is simpler due to the superfunction 
approach. 

Next, one might also wonder how to invert the operator T, i.e., how to compute the above 
polynomials {ga)a- That inversion is the key for the so-called quasi-interpolation approach, and is 
discussed in detail in |2] and [7] (in the PSI context; our superfunction approach already reduced 
the problem to that setup). At base, we seek a simple linear functional /x such that p,* inverts on 
either the convolution tp* or the map f f *' fi. 

Among the various methods, we describe a general recursive approach (see |21; II^)- 
this end, we need first to present this approach in the nondegenerate PSI case, i.e., when the (single) 
generator satisfies the condition '(/^(O) 0. The superfunction method will allow us then to lift 

the result to the FSI setup. 
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Proposition 31 Let 'll: be a compactly supported distribution with ^(0) = 1 that provides approxi¬ 
mation order k in some W 2 , s G IR. Define the polynomials Pa, ot G |a| < k, by the recurrence 

ga-={T(3)g[i, (27) 

I3<a 

where 

c{l) ■={'^ *' ()^)(0) = {fi* ()^)(0) = (O"^ *' V')(0), 7 G TL’l, I 7 I < k. 

Then these polynomials satisfy 

0 “ = V’ * 9a, |a| < k- (28) 

Note that for the expression to rnake sense fi needs to be continuous. 

The other two representations of 0 ( 7 ) are valid for an arbitrary compactly supported distribution 
'Ip- 

Proof. By Theorem El ip*' reproduces all polynomials of degree < k, and hence (cf. e.g., |7j) 
fi *' 0°" = fi * 0" = 0“ *' Ip, Va G 71'^, |q:| < k. 

Thus, c{a) is well-defined. 

Now, given a as above, it is elementary that (since * (f)^ = pj{0) = 1) 

</>»()“= i: w. (rT(o) o'* = 0” + i: cia-mf- 

|/3|<|o| l/3|<|o| 

However, with g^ as in (EH), we obtain (by convolving g^ with assuming by induction that 
'P * 9/3 = 0^ for \(5\ < a, and using the last identity) that 

p * ga--=P * {T - c{a - fi) 0^ = ()“. (29) 

l/3|<|a| 

□ 

Using this proposition with respect to the superfunction p :=Z](/)G# P *' “</>) we obtain the fol¬ 
lowing: 

Theorem 32 Under the assumptions of Corollary the polynomials (pa) from Corollary [23 
satisfy the following recurrence relation: 

9a :=()" - ^ c(a - P) gg, (30) 

P<a 

where 

while 

Here, p * {p is continuous convolution, while a^j, is the finitely supported sequence that appears in 
Corollary (23 Moreover, if each p £ ^ is continuous, we have the alternative discrete convolution 
representation 

c(7,</>) :=(<(>|^^ * ()j^^J(0) = Y PU) apk - j) {-kfi. 

j,kG’Z‘^ 

Remark. Compare the last theorem with Theorem 1 of m- 
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4.5 Good and bad superfunctions 


Every ESI space contains a superfunction. This positive statement can be turned negative: the 
existence of a superfunction in a given ESI space tells us nothing about the structure of the space. 
In contrast, Proposition OHl Corollary EHl and Theorem EH show that the existence of “good” 
superfunctions does lead us to useful conclusions about the space and about the given generating 
set. A particular useful condition is that the Fourier transform of the superfunction be bounded 
away from zero near the origin. This condition is important also from the numerical stability point 
of view. In view of the above, we say that a superfunction ■0 is nondegenerate, if near 0 is 
bounded away from zero. 

We note that the nondegeneracy itself falls short of classifying “good” superfunctions. For 
example, it can be checked that the L 2 -projection of the sinc-function on an ESI space S C L 2 
is always nondegenerate as long as S provides a positive approximation order (by Corollary 
2.6]). However, the superfunctions obtained in this way may prove to be of little use due to their 
slow decay at 00 . Thus, we require a complementary property of a superfunction: we say that the 
superfunction has the finite span property if it is in the finite span of the shifts of the generating 
set <h. Such superfunctions are compactly supported if d> itself is. It is proved in [3] that every 
local ESI space in L 2 contains a superfunction that satisfies the finite span property (with <I> being 
any compactly supported generating set for the space). 

We call a superfunction good, if it is nondegenerate and finitely spanned by the shifts of 4*. Such 
super functions are needed for constructing quasi-interpolation schemes. Indeed, the requirement 
appearing in Theorem EH is exactly that the superfunction ^|J be good. 

Corollary El shows that any function in ^(^(ITI) is of the form (u*<I>)^, for some 27r-periodic 
vector-valued function v. Theorem adds that the vector v associated with a superfunction 
satisfies 


v*G^^sV 


0(1 • 


with k the approximation order of the ESI space ^^(ITI). The known L 2 -theory of approximation 
orders of ESI spaces offers then a recipe for constructing good superfunctions: first solve the 
equation Gv = around the origin (cf. Theorem 1251) and then approximate u by a trigonometric 
polynomial vector u such that v — u has a zero of order k at the origin. This is possible whenever the 
Gramian G is k times continuously differentiable around the origin and G(0) is invertible (see El 
Theorem 4.2]). Next, v — u = 0{\ ■ |^) implies = 0{\ ■ |^), hence (u*4>)^ is a good 

superfunction. 

Once G( 0 ) is not invertible, the notion of a “good” superfunction becomes more subtle. Are we 
only interested in the existence of a superfunction ^jJ £ such that is “reasonably local” and 
'i/’(0) 0, or do we also insist on simple ways to obtain that function from the given generating 

set 4>? Our discussion and development focuses on the latter approach: after all, the SI space is 
given to us in terms of the generating set 4>, and we would like then the analysis to stay as close as 
possible to this set. Once we agree on that principle, it should be clear that “very bad” generating 
sets are not going to yield good superfunctions: for example if 4> is compactly supported and 
<I>(0) = 0, there is no hope to get from 4> in a simple way a superfunction if) with V’(O) 7 ^ 0. The 
ultimate question is how to define “good” vectors <I>. Our suggestion is simple: these are the vectors 
that yield good superfunctions! 


Our next results (in the next subsection) offer analysis of vectors 4> whose Gramian is singular. 
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We show the utility of this analysis by providing a new proof to a famons example of de Boor and 
Hollig concerning the approximation order of C^-cubics on three direction mesh. We then provide 
an example of a “seemingly good” vector < 1 > that cannot yield good superfnnctions. 


4.6 Estimating approximation orders when Gramians are singular 

Theorem HU enables us, at least in principle, to determine the order of approximation provided by 
a given (stationary or nonstationary) FSI ladder in any Sobolev space. But, as we already saw in 
the L 2 -case, such analysis is hard to carry ont if the Gramian | • is singular at the origin. 

The problem is exacerbated by the fact that the entries of may be hard to compnte. 

Let us examine closely the source of the difficulty. For s = 0, the computation of approximation 
orders depends on estimating ratios of the form 

v*G% 

v*Gv 

around the origin. Without loss, one can assnme that the vector v is normalized pointwise. If 
G is continuons at 0 and invertible there, we can then dismiss the denominator, since it does not 
affect the asymptotic behaviour of the above expression. In contrast, if G is singnlar at the origin, 
the denominator might affect the approximation order. The use of the verb “might” is justified: 
roughly speaking, there is hope that the specific vectors v that minimize the numerator are far 
enough from the kernel of G(0). Whenever this is the case, the problem is reduced to examining 
the behaviour of the numerator only. The cnrrent subsection translates the above discussion into 
rigorous analysis. 

We first provide below a theorem that establishes an npper bound on the approximation order 
of an FSI space. The upper bound does not reqnire the invertibility of the associated Gramian. To 
this end, we denote by /Omin(^) the smallest eigenvalue of a positive-definite Hermitian matrix A. 

Theorem 33 Suppose $ C Loo(II) for some neighborhood 11 of the origin. Given any set I C 
2'k' 7L'^\0, denote by k{^,2,s) the order of the zero that the scalar function 

^ ^ Pm\ni.A{u}y)A{u}') '.= 'y ) -|- q ;) 4 *(w -|- ol)*\u) -\- 

OL^X 

has at the origin. Then the approximation order provided by the FSI spaee ^^(Wl) is no larger 
than s)/2. 


Proof. Snppose that 5$ provides approximation order k in W^. Then the characterization from 
Theorem 1^ implies that, for uq : cu vq{iu) that minimizes (12111 pointwise, the expression 


1 

\2k—2s 


=: Mis,,, 


\2k 


is bounded in a neighborhood fl of the origin. Using the identity 

and the fact that k > s, we conclude that VqG^^sVq is bounded above by a constant multiple of 
I . 12 s (compare with the argument used in the proof of Pronosition 111 jl . Thus, 


inf ■ 


v v*Gis>^sV VQGq,^sVQ 


> const ■ 




|2s 


> const inf 






> const 


I-2s 


inf „v. 

hl=i 
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The last inequality uses the assumption that is bounded around the origin. We conclude then 
that, for some C > 0 and a.e. around the origin, 

V*G% on 

inf^^>C|-|-2V^in(Ga,J. 

^ v*G^^sV 

But s~^ (pointwise) a nonnegative definite Hermitian matrix, hence, pointwise, /?min(G$ *) > 
Pmin(^)- The desired result then follows from Theorem I21L □ 

As alluded to before, we know quite precisely when the above upper bound matches the asso¬ 
ciated approximation order. 

Theorem 34 Suppose 4* C Loo(f^) for some neighborhood n of the origin. Let vq he a normal¬ 
ized eigenvector of G% ^ associated with the minimal eigenvalue (i.e., for a.e. uj ^ Ll, the pair 
^Pmm g{uj)),VQ{u))) is an eigenpair of G%g{Ld)). If |nQ4>| is bounded away from zero almost 
everywhere in Ll, then the approximation order of is exactly /c(4>, 27rZi'^\0, s)/2. 


Proof. One only needs to show that the approximation order of 5# is bounded below by 
A:(4>, 27rZi'^\0, s)/2. But 


inf 


v*Glv 


< 


'^oG^c'fo 


Pmm(G^g) 


Pmin(G^ ) 

< -^— < const 


I-2s 


V V*G^,sV V^G^,sVo |nS4>|2| • |2s + J |^*^|2| . |2s 

Theorem 1^ then yields the requisite lower bound. 


Pmin(0(j, 


4.7 Example: bivariate Gbcubics 

The results of the last section raise two questions. The first is whether the upper bounds provided 
in Theorem EHl are useful, i.e., whether they can be applied to solve a nontrivial problem. We 
provide in the current subsection an affirmative answer to this question. 

The other, more fundamental question, is whether the setup of Theorem 1341 is universal, i.e., 
whether we can always dispense with the denominator in the characterization provided in Theo¬ 
rem EH This question is intimately related to the existence of good superfunctions. In the next 
subsection we provide a (-n unfortunate) negative answer to that second question. 

As said, we describe now an example where Theorem EHl applies in a situation when direct 
evaluation of the approximation order is quite complicated. We choose the notorious example of 
an FSI space that reproduces all polynomials of order < 3, but provides only approximation order 
3. The example first appears in |^. A second, completely different, proof of this result appears in 
[3]. Our proof is thus the third one for this result. 

Consider the following two bivariate compactly supported piecewise polynomial functions whose 
Fourier transforms are given by 


4>iiu,v) = i 


(uvw)'^ 

^2{u,v) = 


where w = u-\-v. These functions are known as the Fredrickson elements. With 4> C the 2- 

vector consisting of the above functions, the Gramian is invertible in a punctured neighborhood 
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of the origin. Hence it is still possible to compute enough coefficients in the Taylor expansion of 
1 — to find the first nonvanishing nonconstant term. This complicated analysis was carried 

out in [^. It shows that the first nonzero term in that expansion if of order 6, so provides 
approximation order 3. 

We use here, instead, Theorem ESI to arrive at the same conclusion more easily. Let 2 = 
{(0, 27r), (27r, 0)}. Then 

$((u,u) + a)$((ti,u) + a)* = ^i{u,v)^l{u,v) + ^ 2 {u,v)^ 2 iu,v) + o((|n|^ + |up)^) (31) 


where 




(1 - e-™)(l - e-™)(l - 
{u + V + 27r)2 


'L2(m, v) 


(1 - e-‘“)(l - e-'"’)(l - 
{u + V + 2ttY 


—27r-\-7riu-\-2niv-\-iuv/2-\-\v^ / 2 

27r—7riu-\-iuv/2+iv‘^/2 

ti(?j+27r)^ 

2n—7riv-\-iuv/2+iu^/2 

v(u-\-27ry^ 

— 27r-\-7riv-\-27riu-\-iuv / 2+m^ / 2 
v{u+27pp 


The trace of the matrix + ^^ 2^2 order 4, whereas its determinant is of order 10, so its 

minimal eigenvalue vanishes to order 6 and its maximal eigenvalue to order 4. Since the (matrix) 
terms that were left out of the computation are all of order o(| ■ |^), the eigenvalues of the left hand 
side of m are also of order 6 and 4, respectively. Now Theorem 1331 implies that the approximation 
order of 5$ is at most 3. 

The fact that the approximation order is at least 3 is trivial: the sum <pi + (1)2 yields a super¬ 
function which is nothing but the box spline M 2 , 2,1 (whose approximation order is indeed 3). The 
vector is thus an example where the singularity of the Gramian does not preclude the existence 
of a good superfunction. 

From our standpoint, the G^-cubic vector <I> is “good”, since the finite span of its shifts contains 
a good superfunction. The notoriety of this case is due to the difficulty in asserting that this 
<1> provides approximation order no higher than 3. The fact that the space reproduces all cubic 
polynomials is a sad, misleading, accident. The reader may claim that we ignore the fact that 
here provides an approximation order which is “disappointing”. While that might be the case, it 
goes beyond the realm of this article: we are only interested in ways to capture the approximation 
order of the given space, and not in the construction of SI spaces that provide “satisfactory” 
approximation order. 


4.8 Good and bad superfunctions, continued 

We will now show an example of a vector <1> whose entries seem to be “reasonable” but nonetheless 
it does not admit a good superfunction. This example of a bad together with the example from 
the last section of a good <I>, illustrates the depth of the difficulty in pinning down the notion of a 
good generating set for an FSI space. 

Let <7 be a compactly supported bivariate function whose Fourier transform g has a zero of 
order A: > 2 at each of the 27rK‘^\0 points. Moreover, we assume that 1 — g has a zero of order k 
at the origin. There are many ways to construct such a function. For example, one can take the 
univariate B-spline of order k, apply a suitable differential operator p{D) or order k — 1 to it, and 
then use its tensor product in 2 dimensions. In this case, g is piecewise polynomial of local degree 
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/c — 1 in each of its variables, and with support [0,/c]^. It provides approximation order k (in L 2 , 
for example). 

Now, let e be the bivariate exponential with frequency (27r,0), i.e., e : x e-> We define 

a vector <I> with two components 

+ (j )2 := g - g. 


Here, to recall, D°‘ is the normalized monomial differentiation, viz., 2D^‘^’^'^g is the second derivative 
of g in the first argument. Despite the fact that each provides only approximation order 2, we 
contend that the FSI space 5$(L2) provides approximation order k. We construct, to this end, a 
compactly supported superfunction as follows. 

We choose a vector-valued function v with two components that are trigonometric polynomials 


such that 



around the origin. We then note that the Taylor expansion of order k of ^ around the point (27r, 0) 


_() 0,2 

() 2,0 


At the same time, the Taylor expansion of order A: of ^ around any point of 27rK'^\0 other than 
(27r, 0) is zero. From this, we conclude that the compactly supported function xjj that is defined by 


has a zero of order k + 2 at every point of 27rK‘^\0. Finally, at the origin, ^ — | • \‘^g has a zero of 
order k + 2. 

In order to determine the approximation order provided hy 'ip, we consider the expression 

^ ['4’,$] - I • \^\d? 

|^|2 l-llffp |^|2 

The term 

I • lyp 

|^|2 

is bounded around the origin. The other term, 

[4p,ip] - |^|2 

I • 

remains bounded (around the origin) even when multiplied by | • Thus, V’ provides approxi¬ 
mation order k (in L 2 ), a fortiori <I> provides that approximation order. 

Note that the superfunction xp does not satisfy the desired condition '0(0) ^ 0. In fact, this is 
necessary in a certain sense. Indeed, let r be a 27r-periodic vector-valued function that is continuous 
at the origin and does not vanish there. Let us further define a function / by 


/ = r*$. 
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Then, up to a non zero multiplicative constant, the low order derivatives of / at (27r, 0) coincide 
with the derivatives at the origin of the function 

Since we assume r not to vanish at the origin, it is clear that some second order derivative of the 
above expression does not vanish at the origin. As such, / cannot provide approximation order 
larger than 2. 

While the vector in this example does not yield a good superfunction, it satisfies the following 
positive property: we could use the truncated Gramian g is order to determine the approxima¬ 
tion order of 5$. Indeed, if we normalize the given vector v (i.e., redefine it pointwise as u/|u|), we 
obtain a vector for which yields the correct decay rate {k) at the origin. This means, in 

turn, that the smallest eigenvalue of g still determines the approximation order of the space . 
The superfunction that we obtain in this way (i.e., by using the normalized v) is still not good: it 
decays painfully slowly at oo. 

We close this section with two comments: 

• We do not know at present of an example where the smallest eigenvalue of the truncated 
Gramian G^ ^ does not determine the approximation order of the space provided, of course, 
that $(0) 7 ^ 0. 

• The above example (i.e., of a case when the smallest eigenvalue of G^ ^ determines the 
approximation order while there exists no good superfunction) is very much a multivariate 
phenomenon. It is not hard to prove that such a case is impossible in one variable, and we 
leave it as an exercise to the interested reader. 

4.9 An application: approximation orders of smooth refinable functions 

We provide in this section one of the most interesting applications of superfunction theory: lower 
bounds on approximation orders of smooth refinable vectors. We note that approximation orders 
of refinable vectors are treated in more detail in the next section. However, the current topic fits 
better into the realm of this section. 

At base, our result will show that once is refinable, and once 5$ contains a single nonzero 
function ^|J from a certain class, the stationary ladder generated by must provide an approximation 
order that corresponds to the class of V'- Our definition of the “class” in question requires the Fourier 
transform of V' to decay (in a weak sense) at a certain rate. 

This problem has rich history in the context of PSI ladders (see the introduction to ED)- A 
substantial treatment of the FSI case is given in m- However, that treatment is carried out under 
the assumption that the Gramian G$ is invertible at the origin. In contrast, we focus in this paper 
on the situation where there are multiple solutions to a single refinement equation, and in such 
a case the Gramian of any particular solution is not invertible at the origin. This understanding 
was our motivation to look for an alternative approach to that of m- It is useful to stress that, 
in general, refinable vectors that contain a smooth (even analytic!) function need not provide any 
positive approximation order at all. (An example of this type can be found in (HJ-) Thus, one 
must impose certain side conditions either on the vector $ or on the function ip. 
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Let P be an r X r matrix whose entries are 27r-periodic and measurable. Let be a vector-valued 
function with r components whose entries are in W 2 for some s € IR. We say that <h is refinable 
if the functional equation 

$(2-) = P$ (32) 

is satished. 

Our goal is to prove the following result. In the result, as well as elsewhere in this subsection, 
we use the following notation: 

A :={co € IR'^ ; < |w| < p}. (33) 

Here, p G (0, tt) is arbitrary, but fixed. 

Theorem 35 Let s < 0, and let <I> C W 2 he a solution to HE) . With A as in HE) . assume that 
there exists f G S^{W 2 ) with the following properties: 

1. I/I is bounded above as well as away from zero on A. 

2. The numbers 

\f{-+ o:)?\\l^{A), me 71+ 

QG2"*(27rS'^\0) 

satisfy Xm = 0(2“^™’^), for some positive k. 

Then S'$(hL|) provides approximation order k. 

We approach this result via the notion of the dual equation 

v*{2-)P = v*. (34) 

Here, u is a vector-valued function with r components. We require the dual equation to be valid 
in some (small) ball U centered at the origin, and define the entries of the dual vector v to be 
equal to 0 outside [—vr, 7r]'^\C/. We then extend u to a 27r-periodic vector. Thus, v is supported on 
U + 2^77^, and satisfies there. 

We collect in the next lemma a few simple facts about dual vectors. 

Lemma 36 Let A be as in /T^) . Given any vq defined on A, equation ean he solved on the 
punctured disk 

U :={uj : 0 < Iwl < p} (35) 

so that the solution v satisfies v\a = uq. Moreover, we have then, a.e. on U, 

v*(uj/2^)^(lo/ 2^+ a) = v*(uj)^(uj + 2^a), allmGffi+, a G 27rK'^. 

Proof. We dehne v by v*{ui) :=v*{2u!)P{uj), for all 00 G 2^A, I = -1,-2,_ Then v clearly 

satishes (EU (on U, and hence on U + 2 tt77'^). 

The second part of the lemma is obtained by iterating m times with 

z;*8(. + a) = v*{2-)P^(- + a) = v*{2-)P{- + a)8(- + a) = 

u*(2-)$(2(- + a)). 


□ 
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Proof of Theorem 1231 By Corollary 121 / = for some 27r-periodic r. Denoting by no the 
restriction of r to A, we extend no to a dual vector n by Lemma 1221 Defining '0 by 

•0:=n*8, 

we have, by the same lemma, that, for a.e. on A, and for every nonnegative integer m, 

|^(u;/2-)| = |/M|. 

Thus, in view of our assumptions on /, we conclude that {i/jI is bounded between two positive 
constants around the origin. 

Next, we prove that S^{L 2 ) provides approximation order k. The argument will show, as a 
by-product, that V' G T 2 . 

Since V’ is bounded away from 0 around the origin, it remains to prove, in view of Result |21 
that [^,^]°| • 1“^^ is bounded around the origin (with :=[-0,-0] — |'0|^). Let a; G A, and m a 

positive integer. Then, by the definition of tp and Lemma 1221 

[^,^]°(a;/2™) = ^ + a)= ^ + 2^a) < \m < < C\oj. 

aG27rffi‘*\0 oe27rffi‘^\0 

Thus, = 0(1 • on the punctured ball U of radius p centered at the origin. Since ^ is 

supported on U + 27rK'^, it follows that ip, hence ip, lies in L 2 . Result |21 then applies to show that 
Ip provides approximation order A: in ^ 2 - 

On the other hand, ip = with v measurable and 27r-periodic. Since ip G L 2 , and s < 0, we 
have that ip G VLI- Corollary |21 shows then that ip G ^'^(ITI)- Now, Proposition 1111 implies that ip 
provides approximation order k in W 2 , hence provides also approximation order k in that space. 
□ 

Discussion 37 As the proof of the theorem shows, there is in faet more freedom in the choice of 
A. It suffices to assume that A is compaet, that the intersection ^n2^ has measure zero, and that 
the union contains a neighborhood of the origin. The proof remains essentially the 

same. 

5 Vector refinement equations 

In our studies so far, we considered SI spaces one at a time. There are situations, however, where 
several different SI spaces may stem from one common source. In cases of this type, it is important 
to study the resulting SI spaces in a cohesive, combined, way. 

The best examples of this type are the multiple vector-valued solutions to refinement equations, 
and this is, indeed, the topic of the current section. Let us start with the requisite definitions. 

Let P be an r X r square matrix whose entries are, 27r-periodic (measurable) functions (defined 
on IR'^). The functional equation 

$(2-) = (36) 

is a vector refinement equation, P is a refinement (matrix) mask, and a solution is a 
refinable vector. Here, the entries of the vector <I> are (measurable) functions, or, more generally, 
tempered distributions, defined on IR'^. The rows and columns of the matrix P are, thus, indexed by 
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either the integers 1,..., r, or, more directly, by the entries of ‘h. In this generality, the equation dsni) 
has, as a rule, inhnitely many linearly independent solutions. Indeed, if P is nonsingular around 
zero, then a solution <I> can be chosen arbitrarily on a set A of the ‘dyadic annulus’ type introduced 
in Discussion EZI and then continued to the rest of the Fourier domain using the recipe 

$( 2 u;):=P(a;)$(a;), a; G 2^.4, j = 0 , 1 , 2 ,... 

$(a;) :=P“^(a;)$( 2 a;), uj€ 2^A, j =-1,-2, _ 

Most of the solutions of the above type will decay very slowly (will not be even in Li(lR‘’*)). In 
contrast, if we assume the entries of P to consist of trigonometric polynomials, and if we corre¬ 
spondingly insist on compactly supported solutions, then the solution space is necessarily finite¬ 
dimensional (as explained in detail in the next section). The special instance when the compactly 
supported solution space is one-dimensional is quite well understood (see, e.g., 0,111, [m, and ESI). 
We are therefore primarily interested in the case when there are multiple (in a nontrivial sense) 
compactly supported solutions to the equation PH|). We denote by 

R[P) 

the linear space of all the solutions of (|36|) whose entries are compactly supported distributions. 

The core of study here is the connection between properties of the refinement mask P and its 
corresponding solution(s) <1>. Shift-invariant spaces enter the discussion in a very natural way. For 
example, if the solution vector < 1 > lies in L 2 , then one has the inclusions 

5$(L2) C V{S^{L2)), 

with V the dilation operator / 1 —> /(2-). Due to the above inclusion, we refer to the SI spaces 5$ 
generated by a refinable <I> as a refinable SI space. 

We start our study in this section with the problem of existence of compactly supported solutions 
to (Ihtil) . Our second, and main, topic is the characterization of the approximation orders of the FSI 
space generated by the solutions d> to the refinement equation. This study is based on the premise 
that, in the case where multiple solutions to the same equation exist, the objective should be the 
interplay among those solutions, and not only the individual properties of each one of them. In this 
course of study, we introduce the notions of the combined Gramian and the coherent approximation 
order and connect them with the (i) the approximation orders of the SI spaces generated by the 
solutions to the equations, (ii) the polynomial reproduction property of the mask P, and (iii) the 
sum rules satisfied by P. Finally, we already provided (in S 21 ) lower bounds on the approximation 
order of a refinable SI space in terms of the smoothness of the smoothest function in that space. 

5.1 Compactly supported solutions to the refinement equation 

The structure of the compactly supported solutions of (ESI) was first completely described in ESI- 
We now restate the main result of that paper and provide a different proof for it. We use the partial 
order < on defined by 

a < b a — b ^ TL^^. 

Also, given a nonnegative integer N, we set 

:='[q; G : |(r| := -(-... -1- ^ 

Finally, we recall that the definition of the monomial differential operator D'^ includes the normal¬ 
ization factor 1 / 7 !. 
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Theorem 38 Given an r x r-matrix P whose entries are trigonometric polynomials, set 


N :=max{n : 2”' G spec(P(0))}. 


Then the map 

is a bijection between the collection R{P) of all compactly supported solutions of and the kernel 
ker L of the map 

L-.€^ xZn^€^ xZM-.{w^)a&Z^^{2\^^W^- ^ {D'^-^PmWg-. a^ZM). 

0</3<a 

Proof. Let be a compactly supported distributional solution to (15^ . and denote Wa :=(L>“$)(0), 
a G Since the vector-valued function $ is entire, the vectors Wa are all well-defined. More¬ 
over, one easily concludes from the relation (I3HB (by applying D°‘ to both sides of that identity, 
expanding the right hand side with the aid of Leibniz’ formula, and evaluating the result at 0) that 
the sequence {'Wa)a^’zi’l solves the infinite triangular system 

2l“ln;„= ^ {D»-^P){0)wp, a^TL%. (37) 

0<(i<a 

In particular, {wa)aeZ;^ & ker L. 

Conversely, let w '■={wa)aeZj,r ^ ker L. Then w extends uniquely to a solution to (|T7|) (in order 
to solve uniquely for Wa in (EH) one needs the matrix 2l“l/ — P{0) to be invertible, which is indeed 
the case for every |a| > N, by our assumption of specP(O)). 

Let II • II be any vector norm on C'’. The operator norm on subordinate to || • || will be 

denoted in the same way. We claim that for some constant ^4 > 0, 

\\wa\\ < all aGZ^. (38) 

Let us see first that (EHl yields the existence of a suitable solution to (EH). 

With (1381) in hand, we define (with ()“ the normalized monomial) 
observe that (each of the entries of) g is entire of exponential type, i.e., it satisfies 

lb(w)|| < all uje€^, (39) 

where | • | denotes an arbitrary norm on We need further to show that each of the entries of g is 
the Fourier transform of a compactly supported distribution, which, by the Paley-Wiener-Schwartz 
Theorem (Theorem 19.3 from |34[ p. 375] and Exercise 7.4 from |3()[ p. 216]) amounts to showing 
that (in addition to El) the restriction of g to IR'^ has slow growth at oo. In order to prove the 
requisite slow growth, we follow an argument from [23] : Denoting 

Cl := sup ||P(0II and C 2 := sup \\g{f)\\, 

we pick uj G IR'^, such that 1 < |a;| < 2. By the construction of g, g{2-) = Pg, and hence, for every 
positive n, g{2'^uj) = P{2'^~^ui) ■ ■ ■ P{uj)g{uj). Consequently, 

||5(2’^a;)|| < CsCr < 
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a bound that evidently establishes the sought-for slow growth. 

It remains to prove (EEl). To that end, we pick Nq G IN so that 

/ Q \ 1 

2-~”l|P(0)l| < 1. (j) < 1. 

Since P is a matrix of trigonometric polynomials, there exists ^4 > 0 such that ||(P"P)(0)|| < 
(A/2)l"l/a!, all a. Moreover, by modifying A if need be, we may assume that A satishes the 
estimate for every a G 

In order to prove PS|) for |a| > A'"o, we may assume, by induction, that PS|) holds for all /? with 
\f3\ < |a|. Then, by (jSZI), 


4|«-/3| a\P\ 


= giH ^ 


< 




l/3|<l«l 


2\^-P\{a- py.pi - a! 


I“l 


hpnrp 

p\|a| ^|«t ^ ^|a| 

- UJ a!(l -2-H||P(0)||) “ 

This proves (1551) . and the proof is thus complete. □ 

The theorem can be extended to rehnement equations more general than (|5U1) . For example, we 
can replace the dilation by 2 by a dilation by any matrix M which is expansive, i.e., its spectrum 
lies outside the closed unit disc. 


Theorem 39 Given anr x r-matrix P whose entries are trigonometric polynomials and an expan¬ 
sive d X d matrix M, set 

N := max{n : 0 G spec{M” — P(0)}}. 


Then, the map 

is a bijection between the collection of all (compactly supported) solutions of the refinement equation 
‘h(M-) = on the one hand, and the kernel ker L of the map 


L : X ^ C'-X Zat : ^ ^ P{0)wp : a e Zn), 

O<0<a 


on the other hand. 


Proof. Analogous to that of Theorem 1381 □ 

If the refinement equation is inhomogeneous, viz., a given function is added to its right-hand 
side, then any solution to it is a sum of its specific solution and a solution to the correspond¬ 
ing homogeneous refinement equation. This allows for generalizations of Theorems |3H1 and oni to 
inhomogeneous equations as well. For the exact statement, see m- 

If P(0) is “regular” in the sense that 1 is the largest dyadic eigenvalue of it, the characterization 
of Theorem EHl is much simpler: every right 1-eigenvector of P(0) gives rise to a solution <I> G 
R{P). However, it is easy to generate examples when the largest dyadic eigenvalue of P(0) is 
> 1: for example one can replace P by 2P. In that case, the solutions in R{2P) are obtained 
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by differentiating suitably the solutions in R{P). We close this section with two results related to 
the current discussion. In the first, we describe a general set up in which the solution space R{P) 
is decomposed into the sum of derivatives of solutions to “regular” refinement equations. In the 
second result, we provide an example when such decomposition does not hold. Since the discussion 
here is somewhat tangential to our main study of approximation orders, we skip the proof of the 
following theorem. 


Theorem 40 Given anrxr refinement mask P, let R{P) denote the space of compactly supported 
solutions to and let N be the maximal integer n for which 2" e spec(P(0)). Suppose, that 

that we can find two r x r function-valued matrices T and P such that: (i) T is analytic and 
invertible around the origin, (ii) the entries of P are trigonometric polynomials, (Hi) the matrix 
T{2-)P — PT has a zero of order N + 1 at the origin, and (iv) the Taylor expansion of degree N 
of P around the origin is block-diagonal and the spectrum of each block evaluated at zero intersects 
the set {2^ : j = 0,..., A^} at no more than one point. Let be in R{P), and assume that each of 
the entries of d> has a zero of order I at the origin. Then admits a representation 

N 

= (40) 

j=l 

with G R{P/2^), 4)j(0) 7 ^ 0, and pj is a homogeneous polynomial of degree j, j = I,..., N. 

As mentioned before, the solution space R{P) does not always have such structure, as the 
following counterexample demonstrates. 


Example 41 For some masks P, the decomposition from Theorem W\ is not valid. 


Proof. Let d = 2 and let the mask P satisfy the following conditions: 



' 1 0 o' 


'000' 

P(0) = 

0 2 0 

0-14 

, (D(°’1)P)(0) = 

0 0 0 

0 0 1 


(/)(LO)p)(o) = 0. 


Let us show that the following inclusion fails: 


^ d. 

span[(D“8(0))|„|<Ar : G R{P),^{0) = 0 ] C ^ span[(L»" 

Here ej is the vector in with 1 in position j and zeros 
equivalent to the fact that (jin|) fails. 

The sequences w in C vat-J (in our case r = 3, d = N = 2) indexed by a, |a| < N, we 
envision as ‘long’ vectors with the components Wa, each of length r, all stacked together in some 
fixed order, e.g., in the graded lexicographic order of the a’s. 

The relation (1411) is, again by Theorem I48L equivalent to the following: 


(41) 

elsewhere. By Theorem EHl this is 


d 

ker Lq C^kerLp (42) 

i=i 
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where 


Lo:€''^(N-i) ^ D^-^P{0)wis, 1 < |a| < TV, 

0<l3<a 


^ I (2l“l/ - Pmw^ - Ee,</3<a D--l^P{Q)wp, if a > e, 

I Wa otherwise, 

Now, (SU fails iff its dual statement 


1 < lal < N. 


ran Lq D ran L* 


(43) 


fails. Here (ESI) is obtained from (US) by taking orthogonal complements on both sides and using 
the property kerH = (ranH*)^, which holds, in particular, for any linear map acting on a finite¬ 
dimensional Hilbert space. 



' 0 ' 


0 

Now let := 

0 

, : = 

1 


2 


-2 


. Then 



’ ''^( 0 . 1 ) ' 


0 


’ ^( 1 , 0 ) ' 


0 


0 


0 

w : = 

0 

— T* 

- ^(OA) 

0 

— T * 

- ^( 1 , 0 ) 

0 


0 


^( 1 , 1 ) 


0 


0 


0 


0 


but, by direct calculation, the vector w is not in the range of Lq 


□ 


5.2 Coherent approximation orders 

The general theory of approximation orders of FSI spaces (§4) focuses on the individual space 
5# and its properties. In contrast, when studying the solutions of the refinement equation (ESI), 
we believe that the focus should be on the interplay among the varions solutions, in other words 
on their “common ground”. An attempt to establish a theory that treats simultaneously all the 
solutions of ESI) should be done with care: it is easy to show that different solutions of the same 
refinement equation may have completely different properties, as the following discussion makes 
clear. 


Discussion 42 For j = l,...,r, let (pj he a (sealar-valued) refinable funetion with (scalar) mask 
Pj. That is, eaeh Pj is 2'K-periodie and 4>j{2-) = pjcpj. Define P :=diag{pi,... ,pr), ^ :={4)jYj^i. 
Then ^ is a refinable vector with mask P. For each fixed j, the vector whose jth entry is (pj 
and all other entries are 0 is refinable with respeet to P. Sinee we may seleet the original refinable 
elements {(pj) in a completely arbitrary manner, it is elear that the different solutions {^j)j to the 
same refinement equation may he very different one from the other. 

This discussion reveals another difficnlty that arises when dealing with different solntions to the 
same refinement equation: with Gj the Gramian of that Gramian is singular at the origin. It 
is well-known that this is not an accident: 
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Result 43 H23|l Let ^ C L 2 be a compactly supported refinable vector with Gramian If 

G$(0) is invertible, then the spectral radius q{P{0)) o/P(0) is equal to 1, 1 is the only eigenvalue 
on the unit circle, and 1 is a simple eigenvalue. 

That is, the Gramian of a refinable function is invertible at zero only if the spectrum of P(0) is 
of a special nature, which, in particular, implies that the refinement equation has a unique solution. 
We note that the analysis of the approximation order of this case (viz., a refinable vector whose 
Gramian is invertible at the origin) is carried out in jnj and and is not among our objectives 
here (although we will recall those results momentarily). 

In order to deal with all the solutions of a fixed refinement equation in a combined fashion, 
we introduce first the notions of the combined Gramian and the coherent approximation order of 
the solutions. Let P be a refinement mask, and let (<hi,..., <1>„) be a basis for the solution space 
P(P) of the underlying refinement equation (I3fil) . Assuming that, for some s € IR and for every 
j = 1,..., n, <l>j C 1 ^ 2 ( 24 '^), we define the combined Gramian Gjp^p^ ,, of the refinement equation 
to be the sum of the individual Gramians: 


Gr{p),s 

i=i 

Although the above definition depends on the particular basis that we choose for the solution space, 
our subsequent analysis of ^ is independent of the basis’ choice for the following reason. Let 

B = (<l>i,..., $„) be a basis for P(P). We consider P as an r x re matrix. Thanks to the identity 

n 

= (44) 

i=i 

we conclude that 

GR(ns= E {BB*\-\^^){- + a). 

A new basis for P(P) can be written as BM, with M an re x re constant matrix. Thus the combined 
Gramian for the new basis has the form 


GRiP),s:= E {BMM*B*\-\^n{- + a). 

ol^27v7L^ 

Therefore, for some constants c, (7 > 0, 


CV*GRf^p)^sV < V*Gr(^p)^sV < Gv*Gr(^p)^sV, 

for any vector v. Using the above inequalities, one can easily check that all our subsequent results 
are independent of the choice of B. We also use the notion of the truncated combined Gramian: 

n 


^R{P), 






Definition. Let P be a refinement mask whose solution space R{P) lies in kUKlR'’*). Let 
Gr(p)^s be the corresponding combined Gramian. We say that P(P) (or, in short, P) provides 
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coherent approximation order k if the following condition holds: there exists a neighborhood 
n of 0 such that 


1 

the function Mp^s,k ■ w mf ——- ’ , , 

\l^\2k-2s y v*Gp(^p)^^{uj)v 


belongs to Loo(f^). 


Here, ^ is the truncated version of the combined Gramian, i.e., the expression obtained after 

subtracting from Gp(^p'j^s the terms indexed by a = 0 (cf. the analogous definition of the truncated 
Gramian from Theorem I21j) . 

While Theorem provides ample motivation for the above definition (specifically, it shows 
that the coherent approximation order coincides with the usual approximation order notion in case 
the solution space of (|36|) is one-dimensional), we note that the coherent notion of approximation 
order does not translate immediately into any clear statement on the approximation order of the 
individual solutions. 


Discussion 44 Let us continue with the example in Discussion \4‘A observe that in the case dis¬ 
cussed there, Gp(^p-^ g = diag(G<^^^s,..., with G^.^g the (scalar) Gramian ofcjij, i.e., [(j)j,4>j] 

in the L 2 -case. It follows easily then that the coherent approximation order matches or exceeds the 
approximation order provided by (j)j (for any value of j). 


In order to advance our discussion, we consider vectors v that realize the coherent approximation 
order k. That is, with H C IR'^ some neighborhood of the origin, 


H 9 u; I—>■ v{lo) € C'" 


is measurable, and, a.e. on fl, 


v*(oj)G^^pjg(cv)v(uj) 

v*(c^)Gp(p)^g(uj)v(cj) 




(45) 


We call such v a universal supervector (of order k). A vector u is a regular universal snper- 
vector if (j45l) can be replaced by the conditions that, near the origin. 


'J*GR[p)^gV 


\2s 


and 


v*G'i 


R(P),s'^ 


V V 


V^V 


= o(i-r: 


A regular universal supervector is clearly a universal supervector. 


(46) 


Discussion 45 The existence of a universal supervector is implied (almost automatically) by the 
definition of coherent approximation order. The proof of this fact parallels the proof of the super¬ 
function existence fTheorem, \2^) and is therefore omitted. The regularity of a universal supervector 
V may he implied by either of the following two stronger assumptions: 

(1) The combined Gramian GR(^p-^^g is invertible a.e. around the origin, and the norm of its 
inverse there satisfies 

\\G-p^g\\ = 0{\-\-^^). 

Indeed, in that case v*Gp,gV > cv*v\ ■ for some positive constant c, since ||Gp^^|| is pro¬ 
portional to the reciprocal of the smallest eigenvalue /Omin(G'_R(p),<i) of ond v*Gp(^p-)^gV > 

PminiGR(^P)^g)v*v. On the other hand, if v is a universal supervector, then 

v*Gp^gV < const u* ^ ^$*1 . < const | • p®, 
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where the last inequality follows from the fact that all ^ ^ B are compaetly supported, so their 
Fourier transforms are bounded around the origin. Therefore, the first, henee all the eonditions in 
are satisfied. 

(2) For one of the solutions $ G > c > 0, a.e. in some neighborhood of the 

origin. Indeed, then 

V* Y, • \‘^"v ~ I • 

and the conditions follow from the fact that v is a universal supervector. 

We now connect among the notions of coherent orders, approximation orders, and regular 
universal supervectors. It is worthwhile to note that the following result does not invoke the fact 
that R{P) comprises the solutions to (IHIH) . We do not even need the fact that the individual vectors 
in R{P) are refinable. 


Theorem 46 Assume that the refinement mask P provides coherent approximation order k in 
W|(1R'^). Then 

(a) Let Sp C W 2 ^(]R'^) be the shift-invariant space generated by R{P) (i.e., it is the smallest closed 
shift-invariant subspaee of that contains each entry of each vector in R{P)). Then 

Sp is an FSI spaee and provides approximation order k. 

(h) Let V be a regular universal supervector of order k that is bounded in a neighborhood of the 
origin. Let ^ he a solution of the refinement equation. Then 

(i) = 0(1 • around the origin. In partieular, the funetion defined by :=v*^ 

satisfies the Strang-Fix conditions of order k. 

(a) If, for some positive c, |n*$| > c a.e. in some neighborhood of the origin, then 5$ 
provides approximation order k. Moreover, with a S^ defined by fj:=v*^, the PSI 
space already provides that approximation order. 


Proof, (a) Let i? be a basis for R{P). Then Sp = Sp, with F any vector that contains all the 
entries from all the vectors b & B. Hence Sp is FSI. 

Now, let V : [—7r,7r]'^ —> be a vector that realizes the coherent approximation order k, i.e., 

a.e. on [—vr, vr]'^. 


v*Gr(p),sV 


(47) 


If follows that for a.e. w G 


[—TT, there exists G H such that 


V* ^s{oj)v{uj) 


(48) 


This allows us to represent [—'K,Trf as the disjoint union of sets 0$, <I> G H such that (EHl) holds 
for every 4> G H and a.e. oo G We need, furthermore, to ensure that these sets are measurable. 
We argue the measurability as follows. First, since v is measurable, so are the functions from the 
left-hand side of Therefore, the function 

/min : W mm ' ’ ' , ' 
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is also measurable. Thus, once we define (r2$) by 




v*{u})G^^s{'^)v{uj) ’ 


G 


we obtain the requisite measurability. 

Now, let T(i>, G i?, be the 27r-periodic extensions of the characteristic functions of 0$, G i?. 
Dehning <l>o via its Fourier transform as follows: 


$0 := X] 


we conclude from Corollary 01 that each of the entries of <ho lie in Sp^s- Consequently, S'$q,s C Sp^s- 
On the other hand, the definition of <ho implies that, a.e. on [—vr, vr]'^, 


v*Gl^,sV 

v*Gq,^^sV 


< c 


2k—2s 


(49) 


This, in view of Theorem on shows that provides approximation order k (in 1 T 2 ( 1 R'^)), a 
fortiori its superspace Sp provides that order. 

(b): The regularity of the supervector v implies that, around the origin, 

= 0(1 . 

v*v ^ 

Since v is assumed bounded, we conclude that 

=0(1 • r), 


and therefore 

v*G%^,v = 0(1 • 

for every G R{P)- This proves the first part of (b)(i), while the second part follows from the fact 
that each of the summands u*| •+ap®(^<h*)(-+a) v (that together make up ^ v) is nonnegative, 
hence has to vanish to order 2k as well. 

As to (b)(ii), the analysis above shows that the function ^|J defined by if :=u*<k satisfies 

[^,^]° = [^,^],-H2|.|2« = o(|.r) 


(near the origin). Since we further assume here that |^| > c > 0 around the origin, we also 
conclude that [if^iffs > (?\ ■ P* there. Thus, Q of Theorem 0 holds, and that theorem implies that 
^^/.(ITI (IR*^)) provides approximation order k. □ 


The first part of Theorem I4til leads to the following conclusion: 

Corollary 47 Let P be a refinement mask and let Sp C 1T|(1R'^) be the corresponding SI space. 
If P provides a coherent approximation order k, then there exists if G Sp for which the PSI space 
S.^ C 1T|(1R'^) provides approximation order k. 

Remark. Note that the combined Gramian Gp^g can be defined for any finite-dimensional 
space R of distributional solutions to the refinement equation Likewise, the notion of (regular) 
universal supervectors makes sense with respect to any such space R . The requirement that R be 
the space of all compactly supported solutions actually plays no role in the results of this section. 
The only condition used is that, for each <1> G i?, its Fourier transform be bounded around the 
origin. Therefore, all results of this section are applicable to this more general setup. 
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5.3 Universal supervectors and sum rules 
5.3.1 Known results: singleton solutions in L 2 (]R'^) 

The characterizations to-date of the approximation power of refinable vectors are confined to the 
L 2 setup, and assume, at a minimum, that the Gramian of the (necessarily unique) compactly 
supported solution is invertible at zero (as well as at several additional points). These characteri¬ 
zations allow one to deduce the approximation order provided by the refinable vector directly from 
the mask. Two relevant notions in this context are Condition and the sum rules. We begin 
with the definition of the former. 

Definition. Given k > 0, we say that the refinement mask P satisfies Condition Zp. if there 
exists a vector v of trigonometric polynomials such that, for each I G E, the vector v*{2-)P — di^v* 
has a zero or order k at ttI, while u(0) 7 ^ 0. Here, 

(50) 


is the set of vertices of the d-unit cube. 

Result 48 (HI; mi ) Let P he a refinement mask, and assume that dim R{P) = 1. Let be 
the unique solution of IM, and assume that 4* C L 2 (]R'^) is compactly supported, and that G# is 
invertible at the origin. Then, for k G IN." 

I. If P satisfies Condition Zp, then 5# provides approximation order k. 

II. If provides approximation order k and if G# is invertible at each point of ttE, then P 

satisfies Condition Zp.. 

We note that the compact support assumption on 4> in the above-quoted result can be weakened: 
the essential needed information is about the behavior of the Gramian around E. We refer to H] 
and for more details. 

Gondition Zp, is written on the Fourier domain. It can be equivalently expressed on the “space” 
domain. The equivalent space-based formulations of Gondition Zp are colloquially known as the 
sum rules. We provide, for completeness, the two frequently used versions of these sum rules. The 
second is taken from (see also M), while the first is borrowed from m 

Result 49 Let ^ C L 2 be compactly supported with trigonometric refinement mask P. Let v be a 
vector of trigonometric polynomials. Then the following conditions are equivalent: 

(a) P satisfies Condition Zp with respect to the current v. 

(b) The pairv,P satisfies the 1st version of sum rules: with {v^) and (Hy) the Fourier coefficients 
of V and P, respectively, 

E E <-pPl+ 2 .Q{l + 27) = 2-'" E v%q{^), leE, qe Up_,. 

(Note that is a vector, P^ is a matrix, and 5 ( 7 ) is a scalar.) 
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(c) The pair v, P satisfies the 2nd version of sum rules: With u" = D°‘v{0), a G 71'^, the Taylor 
coefficients of v at the origin, we have: 

^ l^E, \a\<k. 

f3<a 


Proof. The second version of the sum rules is equivalent to Condition as seen by applying 
D°‘ to v*{2-)P — 6iflV*, expanding the first term by Leibniz’ rule, and evaluating the result at ttI, 
l€E. 

The equivalence of the first version to Condition can be argued as follows: first, recall that 
for a finitely supported s G , its Fourier series s has a zero of order k at the origin iff nfc_i lies 
in the kernel of the functional 

Xs: (s* q){0), 

where s*q is either the semi-discrete convolution or the discrete convolution (the statement is true 
for each of the two choices) of the sequence s and the polynomial q. Let m & E, and let Sm be the 
(vector-valued) Fourier coefficients of the function v*{2-)P{- + nm) — 5m,ov*. Thus, Condition 
tells that, for every q G nfc_i and for every m € E, 

0 = As^(9)= XI 

7,»76ffi'^ 76®'^ 

Fixing some I € E, we can write r] = I' + I + 2a, for suitable a G and V G E. Thus, 

'+l+ 2 aq{‘^'l + 2a + 1' + l)e'^'^^^ — dm ,0 E ^- 7 ^( 7 )• 

I'eE 7W 7 

Multiplying the two sides of the last display by and summing over m, we obtain 

0 = E J2^%Pi'+i+2aqi2^ + 2a + r + l)J 2 - E 

I'eE 7 W meE 7 

Thus, 

0 = 2 ^" ^ v%Pi+2aq{27 + 2a + l)-J 2 V-ffiil)- 

7 ,cr 7 

Replacing 7 by 7 — u finishes the proof. □ 

5.3.2 New results: multiple solutions in LF| 

Our analysis of the multiple solution case is based on drawing a connection between universal 
regular supervectors on the one hand, and Condition Z^ (together with its associated sum rules) 
on the other hand. This approach requires some limited regularity formulated in terms of the 
combined Gramian Gr^s of the space R of refinable distributions (see Theorem 1531 for the precise 
assumption). 

So, let P be an r X r refinement mask and let i? be a finite-dimensional space of solutions to 
lying in some Sobolev space: R C VFI (®''^) (iii the sense that every entry of every vector in R lies 
in 147(1R‘^))- Let A: be a positive number. We consider vector-valued functions v, that together 
with P and R satisfy the following assumptions: 
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Assumptions 50 1. There exists a neighborhood 0 of the origin, such that, a.e. on 0 


||G'i?,s(- + ttOII 


0(1), I G E\0, 

0(1-p*), 1 = 0 . 


( 51 ) 


2. The entries of v are 2TT-periodic and measurable. Moreover, v*v is bounded, and bounded 
away from zero, around the origin. 

3. For some k > 0, and for every G the function 4> defined by 4):=v*^ satisfies condition 

In- 


Note that Assumptions EDI are valid (regardless of the value of k) whenever R contains only 
compactly supported solutions, and u is a vector-valued trigonometric polynomial, provided that 
u(0) 0. 

Theorem 51 Let P be an r x r refinement mask, let R C (IR"^) be a finite-dimensional space 
of solutions to and let v be a vector-valued function, so that P, R and v satisfy Assumptions 
\M with respect to some k > 0. If P and v satisfy Condition Zk, then v*G^gV has a zero of order 
2k at the origin, and, in addition, for every ^ £ R, — (u*^>)(0) = 0(| • \^), provided that 
is smooth at the origin. In particular: 

(a) Ifv*Gji^sV ~ I • P* a-e. around the origin, then v is a universal supervector (with respect to R) 

of order k, and hence R provides coherent approximation order k. This, for example, is the 
case ifGn^s is invertible (a.e.) in a neighborhood of the origin, and ||0^^^|| = 0(| • |“^^) a.e. 
there. 

(b) If, for some ^ £ R, |u*<h| > c > 0 a.e. around the origin, then (a) applies, and we further 

conclude that 5$ provides approximation order k. 

Proof. Let a G 2tt7L'^\0, and let m := m{a) > 1 be the smallest integer for which ajTP 0 2'k'7L'^. 
Let ^ £ R. We prove, by induction on m(a), that u*^>(- -|- a) = 0(| ■ |^). For m = 1, we choose 
I £ E\0 such that 2ttI — a £ AtiTL^. Since 

+ a) = v*P Q + , 

the claim follows from the fact that v*P {2 + vr/) = 0(| • |^). For m > 1, we write 



By Condition Z^, v*P{-/2) — v*{-/2) = 0(| • |^). In addition, since m(a/2) = m — 1, the induction 
hypothesis yields that u*(-/2)‘l>((- + a)/2) = 0{\ ■ |^), too. 

Now, fix <i> G i? and define if := Since v is bounded and £ Wf, we can invoke Corollary^ 
and conclude that £ 5$(1F2^). Since tp satisfies (fTT|) (as stipulated in Assumptions [50ll . and since, 
by our argument above, -^(^ -\- a) = 0{\ ■ |^), for every a £ 27rK'^\0, we see that 

^ W. + a)p|.+op* = 0(|.p‘). (52) 

ae27r®‘*\0 
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However, the left-hand-side in the above equality is ^v, and, hence, by summing H52I) over a 
basis of R we obtain that gV = 0(| • 

Next, with ^ as above, the case / = 0 in Condition together with the boundedness of 
around the origin (the latter is embedded in Assumptions imply that 

^(2-) - ^ = (u*(2-)-P - v*)% = 0(1 • 1^). 

Once we assume ^ to be smooth around the origin, the above implies that ^ — "^(O) = 0(| • |*^), as 
claimed. 

The proofs of (a) and (b) are straightforward, hence are omitted. □ 

Since the case of compactly supported solutions and a trigonometric polynomial v is of central 
importance here, we record separately the statement of Theorem 1511 for this case. 

Corollary 52 Let P be a trigonometric polynomial refinement mask, let R{P) C VHI be the space 
of all compactly supported solutions to the refinement equation and let v be a vector-valued 

trigonometric polynomial. Suppose that P and v satisfy Condition Z^. for some k > 0. Then: 

1. For each G R{P), the (compactly supported) defined by satisfies the Strang-Fix 

conditions of order k, and in addition, = 0{\ ■ \^). Consequently, i/u*(0)^(0) 0, 

then 5#(H2*) provides approximation order k, and G S#( 1 T 2 *) is a corresponding superfunc¬ 
tion. 

2. If V* {0)^{0) 0 for some G R{P), then P provides coherent approximation order k, and v 

is a corresponding universal regular supervector. 

At present, we do not know whether Condition is necessary for the provision of coherent 
approximation order k in case dimi?(P) > 1. The results of this type that we are able to prove 
make strong assumptions on the mask P. Below is one such result. The stringent assumption here 
is that P{0) = I. In what follows, we use the notation ^ introduced in 115.21 for the truncated 

combined Gramian. We prove the result only for s = 0, although it extends to other values of s (at 
a cost of a few technical details and more awkward notation). 

Theorem 53 Let P be an r x r trigonometric polynomial refinement mask and let R{P) be the 
space of compactly supported solutions to IM- Suppose that R{P) C L 2 and that the combined 
Gramian satisfies Assumvtion \5IM. is smooth around each I G E and is boundly invertible 

around each I G E. If P{0) = I, the following conditions are equivalent: 

(a) P satisfies Condition with some vector v satisfying Assumvtions \5I\2 and \5iA3. 

(b) There exists a regular universal supervector v of order k for the space R{P). 

In addition, a regular universal supervector v of order k can be always chosen so that, for every 
$ G RiP), 

(u*$)(0) = 0(1 • 1^). 
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Proof. In view of Theorem EH we only need to prove the implication (b) (a). 

We start the proof by noting the identities 

Gr{p) (2-) = + ttI), (53) 

l&E 

and 

G0p(p)(2.) = PG0p(p)P* + ^ (PGp(p)P*)(- + vrO. (54) 

lGE\0 

The first identity is straightforward (and is quite well-known; cf. [23jb The second one is obtained 
by the subtraction of the identity 

(^$*)( 2 -) = 

from the first one. Here, B is the basis for R{P) that was used to define GjK^py 
Let r be any 27r-periodic vector-valued function that satisfies the condition 

^^R{P)'^ ~ ^(1 ■ 1^^)’ 

near the origin. Then (r*G^^p^r)(2-) = 0{\ ■ p^) near the origin. Thus the evaluation at t( 2-) of 
the quadratic form in the right-hand-side of (El leads to a function which has a zero of order 2k 
at the origin. Since each summand there is nonnegative, it follows that, for every I G P\0, 

P(2.)(PGp(p)P*)(. + 7r/)r(2.) = 0(| • p"). 

However, Gp(^p^ is assumed to be boundly invertible around nl, hence we must have that r*(2-)P(--|- 
ttI) = 0(1 • 1^), near the origin, for every I G P\0. In addition, 

P(2.)(PG0p(p)P*)r(2.)=O(|.|2^). (56) 

Now, let u be a universal supervector. Then, (El is satisfied for t := v, hence v satisfies the 
requirements in Condition with respect to each 1^0. It remains to modify v (if need be) so that 
Condition be satisfied at Z = 0, too. Note that so far we have not used out special assumption 
on P. Still, we already know that (El is satisfied for r := u. 

In order to complete our argument, we assume that P = / -|- 0{\ ■ |^) near the origin. We will 
revisit this condition after completing the main part of the proof. This additional condition, when 
applied to (El leads (once we take into account the boundedness and self-adjointness of Gp^pp to 

P(2.)G0p(p)r(2-)=O(|-|2^). (57) 

Thus, we proved that El implies El, and hence, since El is satisfied for r := u, we conclude 
that 

P(2'^-)Gp(p)u(2”-) = 0(1 • for all n gin. (58) 

Our previous analysis then implies that u*(2"'-)P(- -|- ttI) = 0(| • |^), for every integer n > 1, and 
every I G E\0. 

Now, by forming a suitable finite linear combination of u(2’^-), n = 0,1,..., we can construct a 
vector u such that tt(0) = u(0), while u — tt(0) = 0(| • |^) at the origin. Clearly, u*P(--|-7r/) = 0(|-|^), 
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for every I G E\0. Since both u — u(0) and P — I have a /c-fold zero at the origin, we conclude that 
u*{2-)P — u* = 0(1 • 1^). Thus u satishes Condition Z^. 

We finally contend that there is no loss of generality in the assumption that P — I = 0{\ ■ I*') 
around the origin. Indeed, consider the transformation P i—> T~^{2-)PT, where T is a trigonometric 
polynomial-valued matrix, such that T(0) is invertible. For each oj in some small neighborhood of 
the origin, we have the linear isomorphism defined on {^(w) : ‘h G R{P)} by 

Since each $ G R{P) is entire, the isomorphisms induce a corresponding one between the spaces 
R{P) and R{T~^{2-)PT). Moreover, Condition is invariant under this isomorphism, since the 
vector v*{2-)P — Si^v* vanishes to order k at ttI for each / G FI if and only if so does the vector 
{v*T){2-){T-\2-)PT) - 6i^o{v*T). 

So, if we can show that we can choose a matrix-valued polynomial T (of degree smaller than k) 
so that T(0) is invertible and 

FT = r(2-)+0(1-1^), (59) 

then our claim will follow. To this end, let r(0) = I and let the derivatives (Z?"T)(0) be defined 
inductively, according to the partial order of multi-integers a, as solutions to the equation 

(T>"-^F)(0)(T>^r)(0) = 2l“>(T)"r)(0), 0 < |a| < k. (60) 

0<l3<a 

This system is obtained by differentiating dsni) at the origin and is equivalent to For a fixed 

a, the values {D^T)(0), (5 < a are already chosen, and the coefficient of the term D°‘T{0) is 
2i“l — 1^0 (since F(0) = I). Thus, (IHHl) has a solution {D°‘T){0). Thus, P can be assumed to be 
within 0(1 • I*') of the origin. This completes the proof. □ 

5.3.3 Coherent polynomial reproduction 

We restrict our attention again to the space R{P) of compactly supported solutions to the refinement 
equation (EEl). We show that universal supervectors for R{P) are also ultimately connected with 
polynomial reproduction using the shifts of any compactly supported solution from R{P). In short, 
we show that universal supervectors provide universal polynomial reproduction schemes: 

Theorem 54 Let P be a refinement mask whose space of compactly supported solutions R{P) lies 
in VFI(IR'^). Let v be a vector-valued trigonometric polynomial such that, for some k > 0, any one 
of the following conditions holds: 

1. V satisfies Condition Zk- 

2. V and P satisfy the 1st version of the sum rules 

3. V and P satisfy the 2nd version of the sum rules 
4- V is a regular universal supervector of order k. 
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Let a be the vector-valued sequence of the Fourier coefficients of v*, and let G R{P)- Then, with 
ai,... ,ar the entries of a and 4>i,... ,(j)r the entries of <1?, the map 

r 

^ 4)i *' (oi *' q) =: *' {a *' q) 

i=l 

maps n<fc into itself. The map is surjective (hence degree preserving) if and only if v*{0)^{0) 7 ^ 0. 

Proof. By Result I49L conditions 1 through 3 are equivalent one to the other, and, by The¬ 
orem EH each of them implies 4. Condition 4, in turns, implies that, for any <1> G R(T’), the 
compactly supported function if defined by 'il;:=v*^ satisfies the Strang-Fix conditions of order 
k. The discussion preceding Theorem II71 now implies that the semi-discrete convolution operator 
if*' reproduces polynomials of total degree at most k — 1 and, moreover, preserves the degree if if 
is nondegenerate, i.e., V’(O) 7 ^ 0. But if itself is nothing but the semidiscrete convolution *' v*. 
Since the convolution *' is associative, we see that the function *' [v* *' q) is a polynomial in Ilfc-i 
whenever q G nfc_i; moreover, it has exactly the same degree as q whenever u*(0)<h(0) = if{0) 7 ^ 0. 
□ 
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